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Foreword 
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No.  F  49620-92-J-0423/P1.  The  work  described  in  this  report  was  conducted  in 
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project  manager.  Professors  1.  M.  Daniel  and  J.  D.  Achenbach  are  co-principal 
investigators.  Contributors  to  the  program  include  Drs.  H.  S.  Choi,  H.  J.  Chun, 
C.  L.  Tsai  and  S.  C.  Wooh,  and  Mr.  J.  J.  Luo. 

Northwestern  University 


I.  M.  Daniel 
Professor 
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CHARACTERIZATION  OF  DEFORMATION  AND  DAMAGE  IN 
BRITTLE-MATRIX  COMPOSITE  MATERIALS 

1 .  INTRODUCTION 

The  behavior  of  brittle-matrix  composite  materials  is  intimately  related  to  the 
deformation  and  failure  micromechanisms  including  their  exact  sequence  and 
interaction.  These  in  turn  are  related  to  the  properties  of  the  constituents,  matrix, 
fiber  and  interface/interphase,  as  well  as  the  processing  residual  stresses. 

The  micromechanics  of  stress  transfer  and  fracture  of  brittle  matrix  composites 
has  been  studied  experimentally  and  analytically  by  many  investigators  [1-17]. 
Under  longitudinal  tension  damage  initiates  with  the  development  of  multiple 
matrix  cracking,  since  the  failure  strain  of  the  matrix  is  lower  than  that  of  the  fiber. 
This  is  accompanied  by  partial  debonding  at  the  fiber-matrix  interface  and  fiber 
fractures,  which  strongly  depend  on  the  properties  of  the  interface.  Most  studies  to 
date  deal  with  longitudinal  behavior  of  the  unidirectional  composite.  Very  few 
studies  deal  with  behavior  under  transverse  loading  [10,  11,  18-20].  Although  the 
various  failure  mechanisms  are  known  and  understood  their  relative  magnitudes, 
exact  sequence  and  interaction  and  quantitative  effect  on  the  overall  macroscopic 
behavior  vary  from  case  to  case.  The  various  analytical  models  proposed  to  date 
are  deficient  or  incomplete  or  lack  experimental  substantiation. 

The  objective  of  this  investigation  was  to  study  the  deformation  and  failure  of 
brittle-matrix  composite  materials.  Unidirectional  and  cross-ply  continuous-fiber 
ceramic-matrix  composites  under  uniaxial  tension  were  investigated  experimentally 
and  analytically.  The  various  tasks  of  this  study  are  described  below. 

2 .  EXPERIMENTAL  STUDIES  of  UNIDIRECTIONAL  COMPOSITE 
UNDER  LONGITUDINAL  TENSION  [10-14,  21-24] 

2 . 1  Macroscopic  Characterization 

The  material  investigated  is  .  a  unidirectional  SiC/CAS,  calcium 

aluminosilicate  glass  ceramic  (CAS,  Corning  Glass  Works)  reinforced  with 


continuous  silicon  carbide  fibers  (Nicalon,  Nippon  Carbon  Co.).  A  typical 
transverse  section  of  this  material  is  shown  in  Fig.  1.  Properties  of  the  constituents 
are  tabulated  in  Table  1. 

Longitudinal,  transverse  and  shear  tests  were  conducted  to  characterize  the 
macroscopic  response  of  the  composite  material  and  determine  its  macroscopic 
mechanical  properties.  Typical  stress-strain  curves  to  failure  in  the  longitudinal  and 
transverse  directions  are  shown  in  Fig.  2.  The  Difference  between  the  ultimate 
properties  in  the  two  directions  is  extreme.  Additional  tests  were  conducted  to 
determine  the  coefficients  of  thermal  expansion.  Results  are  tabulated  in  Table  2. 

The  characteristic  features  of  the  longitudinal  stress-strain  curve  are  illustrated  in 
Fig.  3.  The  material  initially  behaves  linearly  elastieally  up  to  the  proportional  limit 
(point  A).  This  is  followed  by  a  nonlinear  region  with  deereasing  stiffness  up  to 
the  "bend  over"  point  B.  Then,  follow  a  region  BC  of  stabilized  stiffness  reduction 
rate  (or  sometimes  a  plateau),  another  region  CD  of  increasing  tangential  stiffness 
(strain  hardening),  and  the  terminal  quasi-linear  region  DE.  These  eharacteristie 
features  are  obviously  related  to  the  microscopic  failure  meehanisms  that  develop 
during  loading. 

2.2  Microscopic  Observations  [10-14] 

Failure  mechanisms  were  studied  by  testing  and  observing  specimens  under 
the  mieroscope  (Fig.  4).  The  damage  development  scenario,  observed  directly  in 
part  and  inferred  in  part,  consists  of  the  following  (Fig.  5): 

1 .  Matrix  craek  initiation 

2 .  Craek  propagation  and  fiber/matrix  debonding 

3 .  Crack  saturation 

4.  Fiber  fractures 

5 .  Ultimate  failure 
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Table  1.  Constituent  Material  Properties 


Property 

CAS  Matrix! 

SiC  Fiber2 

Maximum  Use 

- 

Temperature,  “C  (°F) 

1350  (2460) 

1300  (2370) 

Fiber  Diameter  (jim) 

- 

15 

Density  (g/cm3) 

2.8 

2.6 

CoefBcient  of  Thermal 

Expansion, 

10^/°C  (10^/°F) 

5.0  (2.8) 

3.2  (1.8)  -  4.9  (2.7) 

Elastic  Modulus, 

GPa  (106  psi) 

98  (14.2) 

170  (24.6) 

Tensile  Strength, 

124  (18) 

MPa  (ksi) 

(flexural 

1930 (280) 

1  Larsen  and  Adams  ( 1 989) 

2  Prewo  and  Brennan  (1982),  Daniel  et  al.  (1989) 
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Fig.  2  Stress-Strain  Curves  to  Failure  of  Unidirectional  SiC/CAS  Composite 

under  Longitudinal  and  Transverse  Tension 


Table  2.  Measured  Properties  of  SiC/CAS  Unidirectional  Composite 


Property 

Fiber  Volume  Ratio,  Vf 

Ply  Thickness,  t,  mm  (in.) 

Longitudinal  Modulus,  Ei,  GPa  (Msi) 

Transverse  Modulus,  E2,  GPa  (Msi) 

In-plane  Shear  Modulus,  G12,  GPa  (Msi) 

Out-of-plane  Shear  Modulus,  G23,  GPa  (Msi) 

Major  Poisson’s  Ratio,  V12 

Longimdinal  Tensile  Strength,  Fa,  MPa  (ksi) 

Transverse  Tensile  Strength,  F2t,  MPa  (ksi) 

Longimdinal  Ultimate  Tensile  Strain,  £*iit 

Transverse  Ultimate  Tensile  Strain,  £“2t 

Longimdinal  Coefficient  of  Thermal  Expansion, 
ai,  |i£/°K  (  ii£/°F) 

Transverse  Coefficient  of  Thermal  Expansion, 
a.2,  ii£/“K  (  p£/°F) 


Value 

0.39 

0.38  (0.015) 
127  (18.4) 
112  (16.2) 
52  (7.5)  - 
32  (4.7) 
0.18 
448  (65) 
28(4) 

9  X  10-3 
2.5  X  10-4 

4.05  (2.25) 

4.23  (2.35) 


Out-of-plane  Coefficient  of  Thermal  Expansion, 
a3,  p.£/°K  ( {i£/°F) 


4.32  (2.40) 


CO 

CO 

CD 

cn 


Strain 

A:  Proportional  limit 

B:  Bend-over  point 

BC:  Plateau  or  low  tangential  stiffness  region 

CD:  Region  of  increasing  tangential  stiffness 

DE:  Terminal  reoion 

w 


Fig.  3  Stress-Strain  Curve  of  Unidirectional  SiC/CAS  Composite  under 

Longitudinal  Tension 
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r  Real-Time  Observation  of  Failure  Mechanisms 
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Fig.  5  Damage  Development  in  Unidirectional  SiC/CAS  Composite 

under  Longitudinal  Tension 


Matrix  cracking  was  observed  directly  and  recorded  during  loading  (Fig.  6). 
Matrix  crack  density  is  plotted  versus  average  strain  and  compared  with  the  stress- 
strain  curve  in  Fig.  7.  This  plot  shows  that  matrix  cracking  starts  very  close  to  the 
proportional  limit  and  its  density  varies  nearly  linearly  with  applied  average  strain 
up  to  a  saturation  level.  This  occurs  at  a  strain  of  approximately  0.35  to  0.40%  and 
usually  corresponds  to  a  minimum  crack  spacing  of  seven  to  ten  fiber  diameters. 

Fiber/matrix  debonding  is  very  difficult  to  observe  directly.  Indirect 
evidence  of  debonding,  through  ultrasonic  and  acoustic  emission  measurements, 
will  be  discussed  below.  Whenever  fiber  fractures  were  observed,  they  always 
occurred  at  some  distance  from  the  matrix  cracks  and  most  likely  at  the  end  of  the 
debonded  region  (Fig.  8).  This  is  further  illustrated  by  the  SEM  fractograph  of  a 
failure  surface  (Fig.  9). 

2.3  Ultrasonic  Investigations  (21-23] 

Ultrasonic  techniques  were  developed  for  real-time  monitoring  of  damage 
evolution  in  a  unidirectional  ceramic  matrix  composite  under  longitudinal  tensile 
loading.  A  schematic  diagram  of  the  testing  system  used  is  shown  in  Fig.  10.  A 
pair  of  transducers  is  mounted  in  contact  with  the  specimen  and  the  through- 
transmission  signals  are  monitored  during  loading.  Two  types  of  transducers  were 
used,  longitudinal  (compression)  wave  and  shear  wave  transducers.  Wavespeeds 
and  attenuation  corresponding  to  longitudinal  and  shear  waves  across  the  thickness 
(3-direction)  of  the  specimen  were  determined  and  correlated  with  the  macroscopic 
stress-strain  behavior. 

The  change  in  attenuation  of  the  longitudinal  and  shear  waves,  polarized  in 
the  parallel  to  the  fiber  direction,  as  a  function  of  applied  average  strain  is  shown  in 
Figs.  11  and  12  and  compared  with  the  corresponding  stress-strain  curve. 
Considering  the  previous  microscopic  observations  of  damage  development,  it 
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Multiplication  of  Transverse  Matrix  Cracking  in  Unidirectional  SiC/CAS  Composite  Under  Longitudinal  Tension 
25.7ksi,  e  =  0.15%,  (b)  a  =  28.3ksi,  e  =  0.23%  (c)  a  -  32.3ksi,  t  =  0.31%,  (d)  a  =  57.2ksi,  e  =  0.68%. 
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Fig.  7  Stress  and  Matrix  Crack  Density  versus  Strain  in  Unidirectional 
SiC/CAS  Composite  under  Longitudinal  Tension 


TYPTCAL  FTBKR  BREAK 


Typical  Fiber  Fractures  Near  Matrix  Cracks  with  Fiber-Matrix  Debondin; 


Fig.  9  SEM  Fractograph  of  Failure  Surface  under  Longitudinal 
Tension  Illustrating  Fiber  Breaks  and  Pullout. 
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Fig.  10  Schematic  of  Ultrasonic  System 
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Fig.  1 1  Change  in  Attenuation  of  Longitudinal  Waves  due  to  Induced  Damage 
during  Loading  and  Associated  Stress-Strain  Curve 


-400 


(gp)  ‘ ‘uoiiEnueuv 


o 

o 

CM 


O 

O 

O 


O 

o 

CO 


o 

o 

CD 


o 

o 


o 

o 

CM 


CO 

o 


<D 


CD 

d 


lO 


o  ^ 


CJ 

< 

cc 

CD  H 

d  c/) 


CM 

m 

o 


o 

o 


(!s>i)‘o  ‘ssaais 


17 


Fig.  12  Change  in  Attenuation  of  Shear  Wave  due  to  Induced  Damage  during 
Loading  and  Associated  Stress-Strain  Curve 


seems  that  the  attenuation  and  wavespeed  of  the  longitudinal  wave  are  particularly 
sensitive  to  matrix  cracking  and  relatively  insensitive  to  fiber-matrix  debonding.  On 
the  other  hand,  the  attenuation  and  wavespeed  of  the  transverse  (shear)  wave  seem 
to  be  most  sensitive  to  fiber-matrix  debonding.  All  ultrasonic  measures  (attenuation 
and  wavespeeds  of  longitudinal  and  shear  waves)  seem  to  respond  to  the  last  stage 
of  damage  development  consisting  primarily  of  fiber  fractures. 

The  attenuation  and  wavespeed  of  transverse  shear  waves,  polarized  in  the 
normal  to  the  fiber  direction,  were  also  measured  and  used  to  evaluate  the  extent  of 
fiber-matrix  debonding.  The  varying  transverse  shear  modulus  G23  can  be 
determined  from  the  measured  wavespeed  of  the  shear  wave  during  the  entire 
process  of  loading.  In  addition  to  the  wavespeed  the  attenuation  of  the  wave  was 
measured  and  plotted  in  normalized  form  versus  applied  strain  along  with  the  stress- 
strain  curve  (Fig.  13).  It  is  clear  that  the  normalized  signal  loss  starts  at  the  point  of 
matrix  crack  initiation  and  continuous  at  some  moderate  rate  up  to  the  strain 
hardening  region.  At  that  point,  believed  to  be  the  region  of  extensive  fiber-matrix 
debonding,  the  attenuation  rate  of  the  shear  wave  increases  abruptly. 

2.4  Acoustic  Emission  Study  of  Failure  Mechanisms  [24] 

Acoustic  emission  techniques  were  used  to  monitor  damage  mechanisms 
and  evolution  in  a  unidirectional  SiC/CAS  ceramic  composite  under  monotonic 
longitudinal  tensile  loading.  In  addition  to  the  monotonic  loading, 
loading/unloading  tests  were  performed. 

Three  major  types  of  damage  were  investigated;  transverse  matrix  cracking, 
fiber-matrix  interfacial  debonding,  and  fiber  fractures.  Acoustic  emission  events 
were  monitored  in  real  time  during  loading  and  unloading  and  they  were  correlated 
with  the  corresponding  stress-strain  curves  and  microscopically  recorded  damage 
mechanisms. 
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Fig.  13  Change  in  Normalized  Signal  Loss  due  to  Induced  Damage  during 
Loading  with  Superimposed  Stress-Strain  Curve 


The  test  setup  is  shown  in  Fig.  14.  The  cumulative  acoustic  emission 
counts  are  plotted  versus  strain  and  compared  with  the  stress-strain  curve  in  Fig. 
15.  The  onset  of  AE  corresponds  to  the  proportional  limit  of  the  material.  The 
cumulative  AE  output  increases  linearly  with  applied  strain  up  to  the  strain 
hardening  region,  which  corresponds  to  extensive  fiber-matrix  debonding.  At  that 
point,  the  cumulative  AE  rate  increases  abruptly  for  some  time  and  then  it  drops  to  a 
lower  rate  up  to  ultimate  failure.  The  cumulative  AE  counts  were  also  compared 
with  matrix  crack  density  measurements  (Fig.  16).  The  correspondence  of  the 
various  features  of  the  stress-strain,  crack  density  and  AE  curves  and  their  physical 
significance  are  illustrated  in  Fig.  17. 

The  amplitude  distribution  of  AE  events  and  its  relationship  to  the  failure 
mechanisms  were  investigated.  The  variation  of  cumulative  AE  counts  of  various 
amplitude  ranges  with  applied  strain  is  shown  in  Figs.  18  and  19,  and  histograms 
of  AE  events  of  various  amplitude  ranges  are  shown  in  Fig.  20.  It  can  be  deduced 
that  high  amplitude  signals  correspond  to  matrix  cracking  in  region  ABC  and  level 
off  at  matrix  crack  saturation.  It  can  also  be  deduced  that  the  extensive  debonding 
activity  in  region  CD  produces  medium  to  lower  amplitude  signals. 

Several  loading/unloading  tests  were  conducted  as  shown  in  Fig.  21.  Upon 
reloading,  the  material  follows  the  stress-strain  curve  obtained  under  monotonic 
loading.  This  shows  that  the  loading-unloading  process  does  not  introduce  any 
additional  damage.  Stress  and  cumulative  AE  counts  versus  strain  for  some  loading¬ 
unloading  loops  are  shown  in  Figs.  22-24.  These  records  show  that  AE  increases 
up  to  the  peak  applied  strain  and  no  AE  activity  exists  during  unloading. 
Furthermore,  upon  reloading  during  the  following  loop,  no  significant  AE  output  is 
produced  until  the  previous  peak  strain  has  been  exceeded,  i.e.,  the  material  shows 
the  Kaiser  effect.  This  behavior  continues  for  successive  loading-unloading  loops 
up  to  the  strain  level  corresponding  to  strain  hardening  or  matrix  crack  saturation  or 
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Fig.  14  Setup  for  Acoustic  Emission  Test 
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Fig.  15  Stress  and  Cumulative  AE  Counts  versus  Strain  for  Unidirectional 
SiC/CAS  Composite  under  Longitudinal  Tension 
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Fig.  16  Stress,  Cumulative  AE  Counts,  and  Matrix  Crack  Density  versus  Strain 
for  Unidirectional  SiC/CAS  Composite  under  Longitudinal  Tension 
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A:  Proportional  limit 

AE  initiation 
Crack  initiation 

AB:  Slight  stiffness  reduction  (nonlinear  behavior) 

Constant  rate  AE 

Crack  extension  (continuous  cracks) 

BC:  Plateau  or  low  tangential  stiffness  region 

Crack  multiplication 
Fiber/matrix  debonding 
Continued  constant  AE  rate 
C;  Crack  saturation 

Increased  AE  rate 
(for  low  and  moderate  amplitudes) 

CD:  Region  of  increasing  tangential  stiffness 

High  rate  AE  (low  to  intermediate  amplitude) 
Fiber/matrix  debonding/sliding 
DE:  Terminal  region 

Fiber  fracture  and  pullout 
Reduced  AE  rate 


Fig.  17  Correspondence  of  Features  and  Physical  Significance  of 
Stress-Strain,  Crack  Density  and  AE  Curves 
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Fig.  18  Cumulative  AE  Counts  for  High  Amplitude  Thresholds  versus  Strain 

for  a  Unidirectional  SiC/CAS  Composite  under  Longitudinal  Tension 


Slunoo  3V  ©AiiBinuinQ 


26 


Fig.  20  Histogram  of  AE  Events  for  Various  Amplitude  Ranges 
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Fig.  21  Continuous  Loading  and  Loading-Unloading  Stress-Strain  Curves 
for  Unidirectional  SiC/CAS  Composite  under  Uniaxial  Tension 
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Fig.  22  Stress  and  Cumulative  AE  Counts  versus  Strain  for  Loading  and 
Unloading  of  Unidirectional  SiC/CAS  Composite  under 
Longitudinal  Tension  (Loop  1) 
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Stress-strain  curve 
Cumulative  AE  counts 


Fig.  24  Stress  and  Cumulative  AE  Counts  versus  Strain  for  Loading  and 
Unloading  of  Unidirectional  SiC/CAS  Composite  under 
Longitudinal  Tension  (Loop  5) 


extensive  fiber-matrix  debonding  (region  CD).  When  the  corresponding  strain  level 
(approximately  0.5%)  is  reached  during  a  loading-unloading  loop  (loop  6)  an 
abrupt  change  in  AE  rate  occurs  (Figs.  25,  26).  Furthermore,  upon  unloading 
additional  AE  output  is  recorded  which  is  convincing  proof  of  debonding  and 
frictional  sliding.  The  combined  loading-unloading  stress-strain  loops  and 
corresponding  AE  output  are  shown  in  Fig.  27.  The  above  observations  show  that 
most  of  the  debonding  and  sliding  occur  after  matrix  crack  saturation. 

3 .  ANALYTICAL  STUDIES  OF  UNIDIRECTIONAL  COMPOSITE 
3.1  Shear  Lab  Model  for  Longitudinal  Tension  [12,  13,  14] 

A  modified  shear  lag  analysis  started  earlier  was  further  extended  and 

refined.  It  is  based  on  a  cylindrical  composite  element  containing  a  central  fiber  and 
two  transverse  matrix  cracks.  The  basic  assumptions  of  this  analysis  are  that  (1) 
the  fiber  and  matrix  are  bonded  initially  and  even  after  multiple  matrix  cracking,  (2) 
the  axial  load  carried  by  the  matrix  is  proportional  to  the  difference  in  average  axial 
displacements  of  the  matrix  and  fiber,  (3)  the  variation  of  the  shear  stress  in  the 
radial  direction  is  linear  in  the  fiber  and  inverse  quadratic  in  the  matrix,  and  (4)  the 
influence  of  radial  displacements  in  matrix  and  fiber  are  neglected.  The  stress  level 
at  crack  initiation  is  given  by 

aa  =  l^(Fmt-Cynn) 

where 

Ga  =  applied  stress  in  composite 
Gnn  =  residual  axial  stress  in  matrix 
Fnjt  =  matrix  tensile  strength 
Ejn  =  matrix  modulus 

Ej  =  initial  longitudinal  modulus  of  composite  (Ef  Vf  +  Em  Vm) 

Ef  =  fiber  modulus 

Vf,  Vm  =  fiber  and  matrix  volume  ratios,  respectively 
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Fig.  25  Stress  and  Cumulative  AE  Counts  versus  Strain  for  Loading  and 
Unloading  of  Unidirectional  SiC/CAS  Composite  under 
Longitudinal  Tension  (Loop  6) 
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Fig.  26  Stress  and  Cumulative  AE  Counts  versus  Strain  for  Loading  and 
Unloading  of  Unidirectional  SiC/CAS  Composite  under 
Longitudinal  Tension  (Loop  7) 
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Fig.  27  Stress  and  Cumulative  AE  Counts  versus  Strain  for  Loading  and 
Unloading  of  Unidirectional  SiC/CAS  Composite  under 
Longitudinal  Tension 


Thus,  matrix  cracking  in  the  composite  depends  on  the  matrix  tensile 
strength  and  the  processing  residual  stresses.  The  average  residual  stresses  in  the 
fiber  and  matrix  can  be  determined  from  the  following  equilibrium  and  continuity 
conditions  (Fig.  28): 


AT 

= 

-  1000° c 

Ctm 

= 

5.0  X  10-6/°  C 

ttf 

= 

4.0  X  10-6/°  c 

Em 

= 

98  GPa 

Ef 

= 

170  GPa 

Vm 

= 

0.2 

Vf 

= 

0.2 

= 

52  MPa  (7.5  ksi) 

(2) 

(3) 


The  matrix  strength  Fmt  that  controls  matrix  cracking  in  the  composite  is  the 
"in-situ"  strength  which  is  different  from  the  one  measured  in  the  bulk  material.  The 
in-situ  matrix  strength  is  best  determined  by  recording  the  applied  axial 
(longitudinal)  stress  in  the  composite  at  crack  initiation,  Gci,  and  using  eq.  (1). 

Thus,  for  Ga  =  Cci  =  100  MPa,  arm  =  52  MPa,  Em  =  98  GPa  and  Ei  =  127  GPa, 


F 


mt 


Em 

~  "gj"  ^ci  tJrm 


129  MPa  (18.7  ksi) 


After  crack  initiation,  microcracks  propagate  and  combine  to  form  a  pattern 
of  continuous  transverse  cracks  with  an  average  crack  spacing  I  (Fig.  29).  The 
analytical  model  in  question  gives  the  local  stress  distributions  between  two  matrix 
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Fig.  29  Matrix  Cracking  and  Local  Stress  Distributions  in 
Unidirectional  Brittle  Matrix  Composite  under  Longitudinal  Tension 


cracks  [12  -  14].  Closed  form  expressions  were  given  for  the  distributions  of  the 
axial  stress  in  the  matrix  and  fiber  and  the  interfacial  shear  stress: 


where 

CTinx  =  average  axial  stress  in  matrix 
Ofx  =  average  axial  stress  in  fiber 

Ti  =  interfacial  shear  stress 

X  =  axial  coordinate  measured  from  crack  face 
/  =  crack  spacing 


2  2  El 

a  =  —  •  - 

Atf  EfEmVm 


(6) 


4Gf  Gn 


( 1  -  Vf)  2  ' 


rf(  1  +  Vf) 
1  -  Vf 


Gf,  Gm  =  fiber  and  matrix  shear  moduli,  respectively 
rf,  rm  =  fiber  and  matrix  radii  in  model,  respectively 
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Crack  spacing  prior  to  any  significant  debonding  could  be  obtained  by 

setting 

Cfm  =  Fmt  at  X  =  /  /2  in  eq.  (4); 


2  ,  Fin  +  Fi  Gnn 

I  =—  cosh'^ - 

^  Em  Oa  +  El  ( ^rm  "  Fmt) 


(7) 


The  next  failure  mechanism  accompanying  or  following  matrix  cracking  is 
interfacial  debonding  (Fig.  30).  A  debond  initiation  criterion  can  be  expressed  in 
terms  of  a  critical  strain  energy  at  the  intersection  of  the  matrix  crack  and  the  fiber  or 
a  "nominal"  interfacial  shear  strength  Fis.  The  latter  may  not  correspond  to  an 
actual  shear  stress  because  of  the  singularity  of  the  stress  but  it  can  be  correlated  to 
a  similarly  measured  property,  e.g.,  by  means  of  pullout  or  pushout  tests.  Debond 

initiation  then,  could  be  predicted  by  setting  the  interfacial  shear  stress  ti  at  x  =  0 


in  eq.  (6)  to  Fis: 

^  arf  Vm/EmCta  _  c 

”)  T  - 

where 

x,  (0)  =  maximum  interfacial  shear  stress  (at  x  =  0) 
Fis  =  interfacial  shear  strength 


(8) 


The  corresponding  crack  spacing  at  debond  initiation  is  calculated  as 


/di  =  —  tanh 


A 


a 


Fis 


(9) 
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d  =  debonded  length 


I  =  cmck  spacing 


!  1  1-  !  1  1 


Fig.  30  Matrix  Cracking  and  Debonding  in  Unidirectional 
Composite  under  Longitudinal  Tensile  Loading 
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Stress  distributions  following  debonding  have  been  proposed  for  the 
bonded  and  unbonded  regions  including  the  effects  of  residual  stresses  and 

interfacial  friction  as  illustrated  in  Fig.  31.  In  the  debonded  area(o  <  x  <  d) 


II 

X 

(10) 

^  Ga  2  Xf  X 

fx  -  Vf  ff 

(11) 

Xi  =  Xf 

(12) 

where  Xf  =  interfacial  frictional  stress 


In  the  bonded  area  (d  <  x  <  /  -  d) 


(Em 

\El 


+ 


'2VfXfd 

I  Vjh  rf 


cosh  a{l  12  -  x) 
cosh  a{l  H  -  d) 


(13) 


a,  -(Ef  a  iqJrfEmVmg  2Xfd] cosh a(/ 72 - x) 

afx  -  aa  +  Orfj  +  ^  Ca  arf  j 


(14) 


_arfVn 
2Vf 


Xi  = 


Em^  2Vf Xfd 

E,  V„  r, 


sinh  a{l/2-  x) 
cosh  a{l  12  -  d) 


(15) 


The  extent  of  debonding  can  be  calculated  by  setting 
Xi(d)  =  Fis  (atx  =  d) 


(16) 


An  explicit  relation  for  the  debond  length  d  is  obtained  by  neglecting  friction: 


d  =  i 


a  1.^ 


(17) 


where 


^  = 


_  2Fis  Vf 


Vm  Ejn  Oa  ■*"  El  (Tn 
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The  interfacial  shear  strength  which  controls  debonding  is  difficult  to  determine. 
A  determination  of  this  strength  has  been  proposed  in  relation  to  the  matrix  tensile 
strength  by  assuming  debonding  starts  at  matrix  crack  saturation,  i.e., 

^di  —  ^min 


liL  =  ^  coth 


mt 


Vf 


(18) 


For  the  available  material  properties  it  is  found  that 
Fis  =  179  MPa  (26.0  ksi) 

which  can  be  compared  with  microindentation  test  results  for  the  same  material 
[25]. 

Fis  =  248  MPa  (36  ksi) 


The  longitudinal  stress-strain  behavior  of  the  material  was  predicted  based 
on  the  shear  lag  model  discussed  and  using  the  material  properties  obtained.  The 
result  is  shown  in  Fig.  32  where  it  is  compared  with  the  measured  stress-strain 
curve  [22].  The  agreement  between  the  predicted  and  measured  curves  is  good,  at 
least  up  to  the  point  of  debond  initiation.  The  deviation  may  be  due  in  part  to  the 
fact  that  fiber  fracture  was  not  included  in  the  analysis. 


3.2  Interfacial  Debonding  and  Transverse  Shear  Modulus 

An  attempt  was  made  to  calculate  the  transverse  shear  modulus  G23 
throughout  the  loading  history  based  on  Eshelby's  equivalent  inclusion  method  and 
utilizing  the  debond  length  result  from  the  shear  lag  model  discussed  before  [26-28, 
23].  For  perfect  bonding  between  fibers  and  matrix  the  transverse  shear  modulus  is 
predicted  as: 

^  Gm(  1  -  Vf)  +  G23f  (3  -  4Vm  +  Vf) 

Gm  Gn,[  1  +(3  -  4Vn,)  Vf]  +  G23f  (3  -  4Vn,)(  1  -  Vf)  (19) 
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Fig.  32  Stress-Strain  Curves  for  Unidirectional  SiC/CAS 
Composite  under  Longitudinal  Tensile  Loading  (1  -  longitudinal,  2-  transverse) 


For  the  case  of  complete  debonding,  the  transverse  shear  modulus  G23  is 
given  by  [28]: 

Gj3  (1-V)>P 

"(1-V^(P+4x(1-v,5)v,  (20) 

where 

9  =  (1  +  Vm)  [1  +  2(1  -  Vn,)  Q]  (1  -  Vff  +  [4R  +  (3  -  2v2  -  2Vn.  Vi2f)  Q]  (1  -  Vf)  Vf 
+  4(1  -Vm)(l  -Vi2fV2if)(2  +  Vf)VfRQ 

(l+Vm)Eif 
4(1  -  Vi2f  V2lf)Eni 

(l+Vm)E2f 
(1  -  Vi2f  V2lf)Eni 

%  =  (1  -  Vf)2  +  -rr^[(l  -  Vn,  Vi2f)  Q  +  4(1  -  Vn,  Vi2f)]  (1  -  Vf)  VfR  +  Eff  Q 

In  the  case  of  partial  debonding  along  the  fiber  direction,  it  is  assumed  that 

the  fiber-matrix  interface  around  the  circumference  of  the  fiber  is  completely 

debonded.  The  transverse  shear  modulus  is  then  calculated  by  linear  interpolation 

of  the  moduli  for  the  cases  of  perfect  bonding  and  complete  debonding,  i.e., 

—  fr  ^  -  2d  I  1  2d 
^23  -  ^23 1 —  +  '“'23~ 

or  (21) 

G23  =  ( 1 .  ^2^  2d 
G23  V  G23/  ^ 

where  the  debond  length  d  is  calculated  from  eq.  (17).  The  predicted  shear 
modulus  G23  is  plotted  versus  applied  strain  and  compared  with  the  value  obtained 
from  wavespeed  measurements  of  shear  waves  in  Fig.  33.  The  agreement  is  as 
reasonable  as  can  be  expected. 
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Fig.  33  Change  in  Transverse  Shear  Modulus  due  to  Debonding  during  Loading 


Conversely,  the  ultrasonically  determined  shear  modulus  G23  can  be  used  in 
conjunction  with  the  relation  given  above  to  calculate  the  debond  length  d  as  a 
function  of  apphed  average  strain. 

3.3  Numerical  Studies  of  the  Effects  of  Interphases  in 
Unidirectional  Composites  [29] 

A  fiber  with  a  reduced  interphase  stiffness  contained  in  a  unidirectionally 
fiber  reinforced  composite  gives  rise  to  higher  stresses  at  its  interphase  as  well  as  at 
the  interphases  of  neighboring  fibers.  For  a  hexagonal  array  composite,  detailed 
numerical  results  were  presented.  Three  cases  were  investigated:  (1)  perfect 
composite,  (2)  perfect  composite  which  contains  a  single  fiber  with  lower 
interphase  stiffness,  (3)  perfect  composite  which  contains  a  single  fiber  with 
interphase  cracks  and  lower  interphase  stiffness.  Micromechanical  stress 
distributions  at  the  interfaces  of  the  dissimilar  fiber  and  its  neighboring  fibers  were 
calculated  by  the  use  of  the  boundary  element  method  for  an  hexagonal  array 
composite  subjected  to  either  transverse  tensile  or  compressive  uniaxial  loading. 

The  effects  of  fiber-matrix  interphases  on  the  micro-  and  macro-mechanical 
properties  of  unidirectionally  fiber-reinforced  composites  were  also  studied  for 
composites  subjected  to  transverse  shear  loading.  The  interphases  were  modeled 
by  the  spring-layer  model  which  accounts  for  continuity  of  tractions,  but  allows 
radial  and  circumferential  displacement  jumps  across  the  interphase  that  are  linearly 
related  to  the  normal  and  tangential  tractions.  Numerical  calculations  in  basic  cells 
were  carried  out  using  the  boundary  element  method.  For  the  isothermal  case  the 
results  display  the  effects  of  variation  of  interphase  stiffness  on  the  fields  of  stress 
and  deformation  and  on  the  effective  elastic  constants  for  square  array  and 
hexagonal  array  composites  with  and  without  interphase  failure. 

The  effects  of  fiber-matrix  interphases  and  fiber  volume  ratio  on  the  micro- 
and  macro-thermomechanical  behavior  of  unidirectionally  fiber-reinforced  square 
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and  hexagonal  array  composites  were  also  studied  for  the  case  of  uniform 
temperature  variation  (cooling).  Numerical  calculations  in  basic  cells  were  carried 
out  using  the  boundary  element  method. 


3.4  Cylinder  Deformation  Model 

A  new  model  was  investigated  to  remedy  the  deficiencies  of  the  shear  lag 

model  (see  Appendix  for  details).  The  same  cyhndrical  element  between  two  matrix 
cracks  was  analyzed  as  before  (Fig.  34).  A  form  of  the  displacements  in  the  fiber 
and  matrix  was  assumed  to  satisfy  axial  symmetry  and  zero  shear  strain  at  the  outer 
cylinder  boundary.  A  set  of  governing  equations  was  obtained  by  applying  the 
principle  of  virtual  work.  These  equations  refer  to  a  unit  cell  between  two  matrix 
cracks  with  a  crack  opening  displacement  (COD)  and  debonded  regions  at  the  ends 
(Fig.  35). 

Considering  the  periodicity  of  matrix  cracks,  the  governing  equations  can  be 
integrated  to  yield  the  following  nonlinear  damage-strain  relations: 


CODnl  _ 

1 

(N 

r  7  1 

Ei_ 

^  2X^ZdCODi/Rfj 

_t>21  *’22_ 

£2  +  Vl2  J 

^  El 

where 

l-VfVm  1-V[V„  (Om  I 

bi2  = - 

^  f(vf-Vm)  fvj  1+Vf)  /Gf  \ 

l-VfVm  2{l-VfVm)'Gm  } 

b22  =  1  (l-f)(l+Vm)(l-2vj  ^  Gn, 

2(l-VfVj  '  Gf 
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cylinder  displacements:  u(z),  w(2),  /c(z) 
cylinder  stresses:  N^,  N2,  S,  M^,  M2 
cylinder  strains:  €2, 7,  qi,  q2 


Fig.  34  Cylindrical  Element  of  Fiber  and  Matrix  and 
Assumed  Form  of  Displacements 
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Fig.  35  Unit  Cell  Between  . 

COD  and  De 
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longitudinal  elastic  modulus  of  undamaged  composite, 
major  elastic  Poisson's  ratio  of  composite. 


El  = 

Vi  = 

ai  =  average  longitudinal  stress  in  composite. 

£i ,  £2  =  average  longitudinal  and  transverse  strains  in  composite,  respectively. 

Gm,  Gf  =  matrix  and  fiber  shear  moduli,  respectively. 

Vm,  Vf  =  matrix  and  fiber  Poisson's  ratios,  respectively. 

f  =  fiber  volume  ratio. 

CODm  =  matrix  crack  opening  displacement. 

CODi  =  interface  crack  opening  displacement. 


Relations  (22)  can  be  used  to  predict  the  transverse  strain  £2  in  terms  of  the 
longitudinal  strain  ij.  Assuming  that  the  fiber  and  matrix  remain  in  contact,  i.e., 
CODi  =  0 


then. 

Oi 

_  b2i 

-  tJi 

£2  +  V12  — 

Ul-TT 

El 

b22 

1  El 

or 

where 

n/ 

=  e2  +  Vi2^  =  nonlinear  component  of  transverse  strain 
El 


(24) 


^  =  nonlinear  component  of  longitudinal  strain 
El 


b2i 

b22 


"Poisson's  ratio"  relating  nonlinear  components  of  strain 


The  nonlinear  components  of  the  strains  are  illustrated  in  Fig.  36. 
Substituting  the  known  properties  of  the  constituents  and  the  composite  and  the 


52 


-0.1  0.0 


0.1  0.2  0.3  0.4  0.5 

Strain,  c,  (%) 


Fig.  36  Stress  Strain  Curves  of  Unidirectional  SiC/CAS  Composite 

under  Longitudinal  Tension 
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measured  ai  -  ii  longitudinal  stress-strain  relationship  in  eq.  (24),  one  obtains  the 
transverse  strain  £2-  This  predicted  transverse  strain  agrees  perfectly  well  with  the 
one  measured  directly  (Fig.  37). 

The  model  development  up  to  this  point  does  not  entail  failure  criteria.  It 
simply  relates  a  given  state  of  damage  (matrix  cracking)  to  the  average 
(macroscopic)  stresses  and  strains  in  the  composite.  It  makes  use  only  of  linear 
elasticity  with  the  appropriate  boundary  and  interface  conditions. 

Future  work  of  this  model  will  introduce  failure  criteria  for  matrix  cracking 
and  fiber/matrix  debonding,  residual  stresses  and  interfacial  friction.  The  hoped  for 
result  will  be  a  comprehensive  model  capable  of  predicting  damage  evolution  and 
the  complete  stress-strain  curves  to  failure  as  a  function  of  physical  and 
experimentally  measurable  parameters. 

4 .  FAILURE  MECHANISMS  AND  DAMAGE  EVOLUTION  IN 
CROSSPLY  LAMINATES  [30]. 

The  problem  of  crossply  laminates  has  been  studied  extensively  for  polymer- 
matrix  composites.  The  primary  failure  mechanism  in  these  composites  is 
transverse  matrix  cracking  in  the  90°  layer  reaching  a  Umiting  crack  density.  This  is 
followed  by  longitudinal  matrix  cracking  in  the  0°  layers,  local  delaminations  at  the 
intersections  of  matrix  cracks  and,  finally,  fiber  fractures  in  the  0°  layer  leading  to 
ultimate  failure,  Daniel  et  al.  [31].  In  the  case  of  brittle-matrix  composites  the  type, 
sequence  and  interaction  of  failure  mechanisms  are  different  and  characteristic  of  the 
constituent  material  properties.  Some  experimental  work  on  crossply  ceramic- 
matrix  composites  has  been  reported  by  Sbaizero  and  Evans  [32].  They  tested 
symmetric  [0/90]s  SiC/LAS  laminates  in  tension  and  bending.  However,  the 
observed  failure  modes  were  not  easily  explained.  Subsequent  work  on  crossply 
laminates  dealt  with  a  fracture  mechanics  model  of  a  delamination  crack  formed  at 
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Fig.  37  Prediction  of  Transverse  Strains  of  Unidirectional  SiC/CAS 
Composite  under  Longitudinal  Tension 
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the  tip  of  a  notch  cutting  through  0°  and  90°  layers,  see  Sbaizero  et  al.  [33]  and 
Charalambides  [34].  Mall  and  Kim  [35]  also  conducted  experimental  investigations 
of  failure  mechanisms  in  SiC/CAS  laminates.  They  observed  that  initial  failure 
consisted  of  matrix  microcracks  occurring  and  growing  in  all  plies  in  a  random 
manner. 

In  the  present  work  a  systematic  investigation  was  conducted  of  the  failure 
mechanisms,  their  growth  and  interactions  in  a  crossply  SiC/CAS  laminate  under 
uniaxial  tension  [30].  It  represents  the  first  systematic  experimental/analytical  study 
of  a  crossply  laminate  with  brittle  matrix.  Previous  experimental  work  mentioned 
before  reported  random  or  other  edge  dominated  failure  mechanisms  with  no  clear 
trends.  The  failure  mechanisms  observed  and  recorded  in  the  present  study  are 
substantially  different  from  those  reported  before  and  show  a  clear  and  systematic 
trend.  The  initial  stages  of  failure  in  each  layer  are  interpreted  on  the  basis  of  the 
failure  micromechanics  of  the  unidirectional  material  under  longitudinal  and 
transverse  loading.  The  interaction  of  failure  mechanisms  in  the  0°  and  90°  layers 
and  the  later  stages  of  damage  are  explained  on  the  basis  of  classical  lamination 
theory  as  well  as  on  damage  models  available  for  crossply  laminates.  The  various 
failure  mechanisms  and  damage  states  were  correlated  with  the  overall  mechanical 
behavior  of  the  laminate.  The  main  contribution  of  this  work  lies  in  the  integration 
of  the  micromechanics  of  brittle  matrix  single  layers  and  the  macromechanics  of  a 
crossply  laminate  in  predicting  and  interpreting  the  type  and  sequence  of  failure 
mechanisms  and  their  interaction. 

Crossply  specimens  were  loaded  under  the  microscope  and  the  various 
failure  mechanisms,  their  sequence,  interaction  and  overall  damage  development 
were  observed  in  real  time  and  recorded. 

The  first  stage  of  damage  development  consists  of  microcracks  in  the  90° 
layer,  which  then  develop  into  transverse  macrocracks  (Fig.  38).  These 


56 


Fig.  38  Development  of  Failure  Mechanisms  in  Transversely  Loaded 

Ceramic  Matrix  Composite 

(a)  Initial  Radial  Cracks  Around  Closely  Packed  Fibers 

(b)  Initial  Radial  Cracks  Around  Isolated  Fibers 

(c)  Interfacial  Cracks,  and 

(d)  Interconnection  of  Radial  and  Interfacial  Cracks  and 
Formation  of  Macrocracks 
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macrocracks  increase  in  density  up  to  a  saturation  limit,  with  a  minimum  crack 
spacing  approximately  equal  to  the  90°  layer  thickness.  Thereafter,  transverse 
matrix  cracks  are  generated  in  the  0°  layer  (Fig.  39).  These  cracks  also  increase  in 
density  up  to  a  minimum  crack  spacing  of  approximately  eight  fiber  diameters.  In 
the  third  stage  of  damage  development  there  is  an  interconnection  of  the  two  sets  of 
cracks.  Several  of  the  denser  0°  layer  cracks  connect  with  each  of  the  90°  layer 
cracks  in  a  delta-like  pattern  (Fig.  40).  This  is  finally  followed  by  delaminations 
and  additional  cracking  in  the  90°  layer  prior  to  ultimate  failure. 

Damage  development  was  correlated  with  the  macroscopic  stress-strain 
response  of  the  unidirectional  material  under  longitudinal  and  transverse  loading 
and  of  the  crossply  laminate.  It  was  noticed  that  first-ply  failure  in  the  90°  layer 
occurs  at  a  higher  layer  stress  than  the  transverse  tensile  strength  of  the 
unidirectional  material.  This  may  be  attributed  to  the  constraining  effects  of  the  0° 
layer  which  tend  to  increase  the  in-situ  transverse  tensile  strength  of  the  90°  layer. 

Following  crack  saturation  in  the  90°  layer,  cracking  initiated  and  increased 
in  the  0°  layers  at  approximately  the  same  strain  levels  where  similar  cracking  was 
observed  in  the  unidirectional  material  under  longitudinal  loading.  However,  the  0 
layer  stresses  at  the  various  stages  of  damage  development,  were  higher  than 
corresponding  stress  levels  in  the  0°  unidirectional  specimen  for  the  same  damage. 
This  was  explained  as  a  stiffening  or  strengthening  effect  of  the  0°  layers  caused  by 
the  attached  90°  layer.  Ultimate  failure  of  the  laminate  clearly  showed  the 
synergistic  effects  of  the  0°  and  90°  layers  of  the  laminate. 

Damage  development  in  the  90°  layer  was  compared  with  analytical 
predictions.  A  model  developed  previously  for  crossply  graphite/epoxy  composites 
was  used  to  predict  cracking  in  the  90°  layer  and  stiffness  degradation,  both  in  the 
90°  layer  and  in  the  overall  laminate.  The  predictions  were  in  very  good  agreement 
with  experimental  observations. 
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Fig.  39  Matrix  Cracking  in  0°-Ply  Following  Saturation  of  Cracking  in  90° -Ply 
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40  Interconnection  of  0°-Ply  and  90°-Ply  Macrocracks 
(Formation  of  "Delta"  Pattern) 
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The  various  failure  mechanisms  and  the  entire  process  of  damage  evolution 
were  correlated  with  the  macroscopic  stress-strain  curve  of  the  crossply  laminate  as 
shown  in  Fig.  41.  The  effects  of  the  various  failure  mechanisms  on  stiffness 
degradation  is  illustrated. 

It  was  shown  that  the  overall  damage  development  is  characterized  by 
different  scales  at  different  stages.  The  relevant  scales  are  fiber  diameter  and  fiber 
spacing  for  failure  initiation  in  the  90°  layer;  layer  thickness  for  damage  saturation  in 
the  90°  layer;  and  fiber  diameter  for  matrix  cracking  in  the  0°  layer.  Finally,  there  is 
a  scale  crossing  phenomenon  in  the  interaction  between  the  damage  in  the  90°  and  0° 
layers. 

5 .  PUBLICATIONS  AND  PRESENTATIONS 

Reprints  of  selected  papers  published  or  presented  during  the  reporting 
period  are  included  in  the  appendix. 
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ABSTRACT :  Acoustic  emission  (AE)  measurements  were  correlated  with  macroscopic 
stress-strain  behavior  and  direct  microscopic  observations  of  damage  mechanisms  in  a 
unidirectional  ceramic  matrix  composite  under  longitudinal  tensile  loading.  Matrix  crack 
initiation  and  the  onset  of  AE  events  correspond  to  the  proportional  limit  of  the  material. 
Both  the  crack  density  and  cumulative  AE  counts  vary  linearly  up  to  the  point  of  crack  sat¬ 
uration  Thereafter  follows  a  “strain  hardening”  region  corresponding  mainly  to  fiber- 
matrix  debonding  and  manifested  by  an  increase  of  low  to  moderate  amplitude  AE  activity 
and  leveling  off  of  high  amplitude  AE  signals.  At  higher  strains,  beyond  the  “strain  harden¬ 
ing”  region,  there  are  two  characteristic  parts  in  the  AE  output,  one  related  to  further 
damage  (possibly  debonding  and  fiber  fracture)  and  the  other  clearly  related  to  interfacial 
frictional  sliding.  Upon  unloading  further  AE  output  related  to  sliding  is  observed.  AE  ac¬ 
tivity  of  low  to  intermediate  amplitude  continues  in  the  last  stages  of  loading  and  is 
believed  to  be  associated  with  fiber  fractures. 

KEY  WORDS :  ceramic-matrix  composites,  failure  mechanisms,  damage,  acoustic  emis¬ 
sion,  cracking,  debonding,  SiC/CAS  composite. 

INTRODUCTION 

The  behavior  of  brittle-matrix  composite  materials  is  intimately  related  to 
the  deformation  and  failure  micromechanisms  including  their  exact  sequence 
and  interaction.  These  in  turn  are  related  to  the  properties  of  the  constituents, 
matrix,  fiber  and  interface/interphase,  as  well  as  the  processing  residual  stresses. 
It  is  generally  agreed  that  the  possible  failure  mechanisms  under  longitudinal  ten¬ 
sile  loading  are  matrix  cracking,  fiber-matrix  debonding  and  sliding,  fiber  frac¬ 
ture  and  pullout  [1-10].  The  relative  magnitude,  the  exact  sequence  and  quantita- 
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tive  effect  of  these  failure  mechanisms  on  the  overall  behavior  of  the  material  vary 
from  case  to  case.  The  randomness  of  fiber  spacing  and  the  statistical  variability 
of  constituent  properties  cause  the  various  failure  mechanisms  to  span  certain 
ranges  of  applied  strain  during  the  loading  process  [11]. 

Matrix  cracking  and  some  fiber  fractures  have  been  observed  directly  by  reflec¬ 
tion  light  microscopy  during  loading  [5-7].  However,  optical  observation  of  fiber 
matrix  debonding  is  very  difficult.  Ultrasonic  techniques  were  developed  for 
monitoring  of  fiber  matrix  debonding  in  conjunction  with  an  analytical  model 
[9].  The  need  exists  to  complement  and  confirm  previous  observations  and  con¬ 
clusions  by  different  techniques. 

The  objective  of  this  study  was  to  apply  acoustic  emission  (AE)  techniques  for 
detection  and  characterization  of  damage  mechanisms  and  damage  development 
in  unidirectional  ceramic  matrix  composites  under  longitudinal  tension.  Acoustic 
emission  results  were  correlated  with  direct  microscopic  observations  and  the 
macroscopic  stress-strain  behavior  of  the  material. 

EXPERIMENTAL  PROCEDURE 

The  material  investigated  was  unidirectional  SiC/CAS,  calcium  aluminosilicate 
glass  ceramic  reinforced  with  silicon  carbide  fibers  (Nicalon).  The  specimen  was 
obtained  from  Coming  Glass  Works  in  the  form  of  a  3.12  mm  (0.123  in.)  thick, 
16-ply  unidirectional  plate.  The  average  fiber  diameter  is  approximately  15  /im, 
and  the  fiber  volume  ratio  of  the  composite  is  0.40.  The  specimens  were  cut  into 
76  X  3.2  mm  (3  x  0.124  in.)  coupons  and  tabbed  with  glass/epoxy  (Scotchply) 
for  testing.  The  surface  of  the  specimen  was  polished  for  microscopic  obser¬ 
vations. 

Figure  1  illustrates  schematically  the  loading  fixture  designed  for  loading  spec- 


Figure  1.  Setup  for  microscopic  observation  of  failure  mechanisms  in  composite  specimens 
under  load. 
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imens  under  a  reflection  microscope  [5-7] .  Loading  was  applied  and  controlled 
by  a  pneumatic  cylinder  and  recorded  by  a  load  cell.  The  specimen  was  in¬ 
strumented  with  strain  gages  so  that  the  stress-strain  behavior  was  continuously 
recorded  while  the  specimen  was  observed  on  the  surface  by  a  CCD  camera  and 
recorded  on  video  tape.  The  matrix  crack  density  at  various  applied  stress  levels 
was  then  easily  determined  from  the  video  records  after  the  test  was  finished  by 
measuring  the  average  spacing  between  matrix  cracks. 

Additional  specimens  were  prepared  and  loaded  in  a  servohydraulic  testing 
machine  (Instron).  These  specimens  were  instrumented  with  strain  gages  to 
record  the  stress-strain  behavior,  while  the  acoustic  emission  output  was  moni¬ 
tored  concurrently  with  a  150  kHz  AE  sensor  attached  to  the  specimen  (Figure 
2).  Acoustic  emission  data  in  the  form  of  cumulative  counts  (hits)  above  a  prede¬ 
termined  threshold  were  recorded  and  processed  by  a  system  incorporating  a  ded¬ 
icated  computer  and  a  transient  record  analyzer  (Locan-AT  with  TRA,  Physical 
.Acoustics  Corp.).  Some  of  the  specimens  were  loaded  up  to  failure  in  a  con¬ 
tinuous  fashion  while  others  were  tested  under  a  number  of  loading-unloading 
cycles  with  increasing  maximum  load. 

RESULTS  AND  DISCUSSION 


Monotonic  Loading 

Results  of  the  macroscopic  stress-strain  behavior,  matrix  crack  development 
and  acoustic  emission  output  were  obtained,  analyzed  and  correlated.  Typical 
photomicrographs  showing  the  development  and  multiplication  of  matrix  crack¬ 
ing  are  shown  in  Figure  3.  The  crack  density  (number  of  cracks  per  unit  length) 
was  measured  from  such  photomicrographs  and,  after  normalizing  it  (by  multi¬ 
plying  by  five  times  the  fiber  diameter),  it  was  plotted  versus  applied  strain  in 
Figure  4.  It  was  found  that  matrix  cracking  occurred  primarily  in  the  0.1-0.4% 
strain  range.  Superimposed  in  Figure  4  is  a  typical  stress-strain  curve  to  failure 
under  longitudinal  tensile  loading.  This  curve  displays  several  characteristic 
features  which  are  related  to  the  failure  mechanisms  and  the  failure  process  as 
will  be  discussed  later.  Figure  5  shows  the  acoustic  emission  output  in  the  form 
of  cumulative  AE  counts  plotted  versus  applied  strain.  Superimposed  for  refer¬ 
ence  is  the  same  stress-strain  curve  of  Figure  4. 

The  various  characteristics  and  measured  parameters  of  the  above  results  were 
correlated  over  five  different  regions  of  the  loading  history  (Figure  6).  In  the  ini¬ 
tial  region,  the  stress-strain  curve  is  linear  up  to  the  proportional  limit  A  (tr  =  97 
MPa,  e  =  0.073%).  In  this  region  no  apparent  AE  activity  was  recorded  and  no 
matrix  cracking  was  observed.  As  the  load  is  increased,  the  stress-strain  curve 
becomes  slightly  nonlinear  because  of  the  stiffness  degradation  resulting  from 
matrix  cracking  (region  AB,  a  =  97-213  MPa,  €  =  0.073-0.20%).  Several  ini¬ 
tial  microcracks  of  typical  length  equal  to  one  to  two  interfiber  distances  were 
observed  in  this  region.  Acoustic  emission  output  started  in  this  region  and  in¬ 
creased  linearly. 

The  characteristic  “bend-over  point”  B,  is  followed  by  a  region  BC  of  reduced 


AE  Study  of  Failure  Mechanisms  in  Ceramic  Matrix  Composites 


1949 


Figure  2.  Schematic  diagram  of  acoustic  emission  monitoring  under  uniaxial  tensile  load¬ 
ing. 
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Figure  4.  Stress  and  matrix  crack  density  versus  strain  in  unidirectional  SiC/CAS  coi 
under  longitudinai  tension. 


Figure  5.  Stress  and  cumulative  AE  counts  versus  strain  for  unidirectional  SiC/CAS  com¬ 
posite  under  longitudinal  tension. 
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Figure  6.  Stress,  cumulative  AE  counts,  and  matrix  crack  density  versus  strain  for  unidirec¬ 
tional  SiC/CAS  composite  under  longitudinal  tension. 


but  stabilized  tangential  stiffness  (a  =  213-252  MPa,  e  =  0.20-0.36%).  In  this 
region  the  matriiT  microcracks  extend  into  macrocracks  and  increase  in  number 
throughout  the  gage  section  of  the  specimen.  Both  crack  density  and  cumulative 
AE  counts  increase  at  a  constant  rate  up  to  a  point  C  where  the  matrix  cracking 
is  saturated  (a  =  252  MPa,  e  =  0.36%).  The  minimum  crack  spacing  at  this 
saturation  point  is  approximately  ten  fiber  diameters. 

The  following  region  CD  (a  —  252-306  MPa,  e  =  0.36-0.52%)  is  one  of  in¬ 
creasing  tangential  stiffness  or  a  “strain  hardening”  region.  Although  matrix 
cracking  appears  to  be  saturated,  increased  AE  output  at  a  higher  rate  appears. 
It  is  believed  that  this  AE  activity  is  related  to  fiber-matrix  debonding  and/or  slid¬ 
ing  occurring  in  this  region  as  indicated  by  ultrasonic  measurements  and  analysis 

The  terminal  region  DE  (a  =  306-433  MPa,  e  =  0.52-0.86%)  is  character¬ 
ized  by  a  gradually  decreasing  stiffness  due  to  fiber  fractures  and  pullout.  AE  ac¬ 
tivity  continues  at  a  reduced  but  stabilized  rate. 

The  amplitude  distribution  of  AE  events  and  its  relationship  to  the  failure 
mechanisms  were  investigated.  The  variation  of  cumulative  AE  counts  of  various 
amplitude  ranges  with  applied  strain  is  shown  in  Figure  7,  and  histograms  of  AE 
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Figure  7.  Cumulative  AE  counts  of  various  amplitude  ranges  for  a  unidirectional  SiC/CAS 
composite  under  longitudinal  tension. 
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Figure  8.  Histograms  of  AE  events  for  various  amplitude  ranges. 


events  of  various  amplitude  ranges  are  shown  in  Figure  8.  It  is  seen  that  matrix 
cracking,  which  is  the  primary  failure  mechanism  in  region  ABC,  Produces  a 
broad  band  of  AE  activity.  However,  the  highest  amplitude  signals  (>75  dB)  oc¬ 
cur  only  in  this  region  and  cease  at  point  C,  the  crack  saturation  level  (Figure  ) . 
Following  crack  saturation,  middle  to  low  amplitude  events  increase  at  increased 
rates  over  the  “strain  hardening”  region  CD.  These  results  suggest  that  fiber- 
matrix  debonding  and  sliding,  the  dominant  failure  mechanism  m  region  CD, 
produce  AE  signals  of  lower  amplitude  but  longer  duration.  In  the  final  region 
DE  the  AE  events  of  low  to  moderate  amplitude  are  significantly  fewer  than  in  re¬ 
gion  CD,  while  high  amplitude  signals  (70-74  dB).  begin  to  appear  again.  Tins 
phenomenon  could  be  explained  by  the  fact  that  fiber  fractures  produce  high 
amplitude  signals  of  narrow  pulse  width  but  few  in  number. 


Loading-Unloading 

In  addition  to  the  continuous  loading  tests,  loading-unloading  tests  were  con¬ 
ducted.  Figure  9  shows  that  the  envelope  of  the  loading-unloading  stress-strain 
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curves  matches  the  monotonic  stress-strain  curve.  This  implies  that  the  loading- 
unloading  process  produces  no  significant  additional  damage  in  the  material. 
Figures  10(a)- 10(d)  show  cumulative  AE  count  curves  along  with  the  corre¬ 
sponding  stress-strain  curves  for  the  various  loading-unloading  loops.  Figure  11 
shows  the  records  for  all  loops  combined.  It  is  noted  that  in  the  region  below 
point  C,  or  the  first  four  loops  of  loading-unloading,  no  AE  output  is  produced 
during  unloading  and  reloading  up  to  approximately  the  previous  peak  loading. 
This  confirms  the  classical  Kaiser  effect  and  also  implies  that  no  significant 
debonding/sliding  exists  before  saturation  of  matrix  cracking.  Any  frictional  slid¬ 
ing  during  unloading  would  produce  additional  AE  counts. 

In  the  “strain  hardening”  region  CD  where  fiber-matrix  debonding  is  believed 
to  be  dominant,  the  AE  count  rate  increases  during  loading,  but  no  activity  is 
recorded  during  unloading  (loops  4  and  5). 

In  regions  beyond  point  D  (loops  6,7  and  final  loading)  a  pronounced  abrupt  re¬ 
duction  in  AE  count  rate  occurs  during  loading.  The  AE  output  continues  in¬ 
creasing  at  a  slightly  reduced  rate  during  the  initial  stage  of  unloading.  This  AE 
activity  during  the  last  stage  of  loading  and  the  initial  stage  of  unloading  (in  loops 
6  and  7)  is  attributed  to  frictional  sliding  in  the  debonded  fiber-matrix  interface. 
The  characteristic  rate  change  in  AE  counts  in  the  loading  curves  of  loop  7  and 
final  loading  step  occurs  at  a  strain  equal  to  the  peak  strain  reached  in  Ae  pre¬ 
vious  loading.  It,  seems  that  the  two  parts  of  AE  output,  the  increasing  rate  part 
and  the  reduced  constant  rate  one,  demonstrate  separate  and  independent  Kaiser 
effects.  The  first  part  of  the  AE  output  is  believed  to  be  associated  with  further 
debonding  and/or  fiber  fractures  whereas  the  second  part  is  primarily  due  to  fric¬ 
tional  sliding. 


Figure  9.  Continuous  loading  and  loading-unloading  stress-strain  curves  for  unidirectional 
SiC/CAS  composite  under  uniaxial  tension. 


Strain,  t,  (%) 


Figure  10.  Stress  and  cumulative  AE  counts  versus  strain  for  loading  and  unloading  of  uni 
directional  SiC/CAS  composite  under  longitudinal  tension. 
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Figure  10  (continued).  Stress  and  cumulative  AE  counts  versus  strain  for  loading  and 
unloading  of  unidirectional  SiC/CAS  composite  under  longitudinal  tension. 
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Figure  10  (continued).  Stress  and  cumulative  AE  counts  versus  strain  for  loading  and 
unloading  of  unidirectional  SiC/CAS  composite  under  longitudinal  tension. 
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Figure  10  (continued).  Stress  and  cumulative  AE  counts  versus  strain  for  loading  and 
unloading  of  unidirectional  SiC/CAS  composite  under  longitudinal  tension. 
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F/gure  1 1.  Stress  and  cumulative  AE  counts  versus  strain  for  ioading  and  unioading  of  uni¬ 
directional  SiC/CAS  composite  under  longitudinal  tension. 


SUMMARY  AND  CONCLUSION 

A  correlation  was  made  between  macroscopic  stress-strain  behavior,  acoustic 
emission  measurements  and  direct  microscopic  observations  of  damage  develop¬ 
ment  m  unidirectional  SiC/CAS  composites.  Dominant  damage  mechanisms  in 
different  loading  stages  were  identified  by  the  correlation. 

Under  monotonic  tensile  loading  matrix  crack  initiation  and  onset  of  AE  events 
correspond  to  the  proportional  limit  of  the  stress-strain  curve.  The  matrix  crack 
density  and  the  cumulative  AE  counts  vary  linearly  with  strain  through  the 
nonlinear  region  between  the  proportinal  limit  and  the  beginning  of  the  “strain 
hardening  region.  At  this  point  matrix  cracking  reaches  a  saturation  level  and 
the  high  amplitude  AE  activity  levels  off  while  the  low  to  intermediate  amplitude 
AE  rate  increases.  The  latter  AE  activity  occurring  over  the  “strain  hardening” 
region  is  believed  to  be  associated  with  fiber-matrix  debonding. 

In  the  case  of  loading-unloading  cycles  it  is  noted  that  the  material  demon¬ 
strates  the  Kaiser  effect  for  loadings  below  the  matrix  crack  saturation  point. 
Beyond  this  point,  increasing  fiber-matrix  debonding  and  fiber  fractures  produce 
changes  in  the  AE  rate  and  deviations  from  the  Kaiser  effect. 

At  higher  strain  levels  beyond  the  “strain  hardening”  region  an  abrupt  reduction 
is  noticed  in  AE  rate  accompanied  by  ftirther  AE  output  during  the  initial  stage 
of  unloading.  This  characteristic  behavior  is  attributed  to  frictional  sliding  during 
loading  and  unloading  in  the  debonded  fiber-matrix  interfaces. 

It  was  shown  in  this  study  that  acoustic  emission  is  a  useful  tool  for  studying 
the  failure  mechanisms  of  composite  materials.  For  complete  identification  of 
failure  mechanisms,  all  three  kinds  of  damage  should  be  observed;  matrix  crack- 


AE  Study  of  Failure  Mechanisms  in  Ceramic  Matrix  Composites 


1961 


ing,  fiber-matrix  debonding,  and  fiber  fracture,  of  which  the  matrix  cracking  and 
fiber  fracture  are  directly  observable  by  the  microscopic  investigation.  In  the 
absence  of  reliable  and  practical  methods  for  direct  observation  of  debonding, 
AE  is  a  good  alternative  method  that  provides  indirect  information  about  fiber- 
matrix  debonding. 
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Real-time  ultrasonic  techniques  were  developed  for  monitoring  damage  in  a  unidirectional  ceramic-matrix  composite 
under  longitudinal  tensile  loading.  Specifically,  shear-wave  transducers  producing  waves  polarized  in  the  transverse  to  the 
fiber  direction  were  used  in  contact  with  the  specimen  to  detect  the  initiation  and  propagation  of  fiber-matrix  debonding, 
and  to  determine  the  transverse  shear  modulus  and  its  degradation.  The  ultrasonically  measured  transverse  shear  modulus 
and  its  degradation  was  in  reasonably  good  agreement  with  a  prediction  based  on  a  modified  shear  lag  model  and 
interpolation  between  available  solutions  for  fully  bonded  and  fully  unbonded  fibers. 


1.  Introduction 

’V 

Ceramic-matrix  composites  are  receiving 
greater  attention  because  of  their  potential  use  at 
high  temperatures.  Composites  consisting  of  a 
ceramic  matrix  reinforced  with  high  strength 
fibers  have  substantially  higher  strength,  higher 
ultimate  strains,  better  fracture  toughness  and 
ductility  than  comparable  monolithic-matrix  ma¬ 
terials.  These  materials  are  attaining  strengths  of 
1400  MPa  (200  ksi)  and  ultimate  strains  of  over 
1%.  Typical  composites  include  silicon-carbide 
matrix  reinforced  with  silicon-carbide  fiber 
(SiC/SiC)  and  glass  ceramics  reinforced  with  sili¬ 
con-carbide  fibers,  such  as  SiC/LAS  (silicon  car¬ 
bide/lithium  aluminosilicate)  and  SiC/CAS 
(silicon  carbide /  calcium  aluminosilicate);  see 
Brennan  and  Prewo  (1982),  and  Prewo  and  Bren¬ 
nan  (1982). 
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ing  and  Applied  Science,  Northwestern  University,  Center  for 
Quality  Engineering  and  Failure  Prevention,  327  Catalysis 
Building,  2137  North  Sheridan  Road,  Evanston,  IL  60208- 
3020,  USA. 


The  macroscopic  behavior  of  the  material  is 
intimately  related  to  the  constituent  properties 
(fiber,  matrix  and  interphase),  the  micromechan¬ 
ics  of  stress  transfer  between  the  constituents, 
and  the  failure  mechanisms  and  their  develop¬ 
ment  under  load.  The  micromechanics  of  stress 
transfer  and  fracture  of  brittle  matrix  composites 
has  been  studied  analytically  by  many  authors; 
e.g.,  Achenbach  and  Zhu  (1989),  Aveston  et  al. 
(1971),  Aveston  and  Kelly  (1973),  Greszczuk 
G969),  Hsueh  (1988),  Lee  and  Daniel  (1992), 
Marshall  et  al.  (1985),  McCartney  (1989),  Stang 
and  Shah  (1986).  Although  the  various  failure 
mechanisms  are  known,  their  relative  magnitude, 
exact  sequence  and  quantitative  effect  on  overall 
behavior  vary  from  case  to  case. 

Failure  mechanisms  have  been  observed  in 
real  time  by  Marshall  and  Evans  (1985)  and 
Daniel  et  al.  (1989a,  1989b,  1992),  but  the  obser¬ 
vations  must  be  made  under  the  microscope  and 
are  limited  to  the  surface  only.  Furthermore, 
some  mechanisms,  such  as  fiber-matrix  debond¬ 
ing,  are  very  difficult  to  observe.  There  is  a  need 
to  bridge  the  microscopic  and  macroscopic  obser¬ 
vations  by  means  of  nondestructive  monitoring  of 
internal  damage.  The  ability  to  monitor  damage 


0167-6636/94/$07.00  ©  1994  -  Elsevier  Science  B.V.  All  rights  reserved 
SSDI  0167-6636(93)E0081-C 


380 


5.C.  Wooh  and  LM.  Daniel  /  Fiber-matrix  debonding 


development  during  loading  without  interrupting 
the  test  is  very  important  because  the  specimen 
must  be  loaded  monotonically  to  failure  at  a 
controlled  loading  rate.  Wooh  and  Daniel  (1992) 
demonstrated  that  ultrasonic  methods  can  be  used 
to  monitor  damage  evolution  in  real  time.  Ultra¬ 
sonic  measurements  of  attenuation  and  velocity 
of  longitudinal  and  transverse  waves  were  mea¬ 
sured  and  correlated  qualitatively  with  failure 
mechanisms. 

Of  the  properties  involved,  the  fiber-matrbc 
bonding  condition  plays  an  important  role  on  the 
mechanical  behavior  of  the  composite.  For  exam¬ 
ple,  if  the  bonding  is  too  weak,  stresses  are  trans¬ 
ferred  poorly.  On  the  other  hand,  if  the  bonding 
is  very  strong,  the  material  behaves  like  a  mono¬ 
lithic  brittle  one.  Fiber-matrix  debonding  under 
stress  is  a  sensitive  indicator  of  material  response 
and  thus  it  is  important  to  investigate  it  nonde- 
structively. 

This  paper  describes  an  ultrasonic  technique 
for  real-time  monitoring  of  fiber-matrix  debond¬ 
ing  and  a  correlation  between  the  ultrasonically 
measured  transverse  shear  modulus  and  an  ana¬ 
lytical  prediction. 


2.  Theoretical  background 

The  material  investigated  was  a  unidirectional 
SiC/CAS  composite,  calcium  aluminosilicate 
glass  ceramic  reinforced  with  silicon  carbide 
fibers.  Figure  1  illustrates  a  typical  longitudinal 
stress-strain  curve  and  corresponding  failure 
mechanisms  for  this  material.  The  curve  shows 
four  primary  regions  of  interest.  The  initial  re¬ 
gion  of  the  curve  (A)  is  the  linear  elastic  range 
where  virtually  no  damage  could  be  found.  This 
linear  region  extends  at  least  up  to  the  level 
where  the  matrix  stress  or  strain  reaches  the 
ultimate  value  of  the  matrix  tensile  strength  or 
fracture  strain.  The  initial  longitudinal  modulus 
of  the  composite  is  obtained  by  the  rule  of  mix¬ 
tures: 

E,=E,V,  +  E^V^,  (1) 

where  Ef,  E^  are  the  moduli  of  the  fiber  and 


Fig.  1.  Schematic  illustration  of  stress-strain  curve  and  failure 
mechanisms  of  unidirectional  SiC/CAS  composite. 


matrix  and  Vf,  are  the  fiber  and  matrbc  vol¬ 
ume  ratios,  respectively.  Assuming  the  ultimate 
matrix  strain  to  be  lower  than  that  of  the  fiber, 
the  stress  level  at  matrix  crack  initiation  is  given 
by  the  relationship 

O’]  =  (2) 

Em 

where  F„„  is  the  matrix  tensile  strength  and 
is  the  axial  residual  stress  in  the  matrbc. 

Following  crack  initiation  and  propagation 
across  the  specimen  section,  accompanied  by 
fiber-matrix  debonding,  the  curve  tends  to  bend 
to  a  plateau  or  reduced  stiffness  region  B  corre¬ 
sponding  to  crack  multiplication  and  further 
fiber-matrix  debonding.  In  the  middle  of  this 
region  (C),  matrix  cracking  becomes  saturated 
while  debonding  continues.  The  curve  enters  re¬ 
gion  D  as  the  interfacial  debonding  is  saturated. 
Fiber  fractures  are  observed  somewhere  in  this 
range  of  stress  and  the  material  finally  fails  glob¬ 
ally. 

Figure  2  illustrates  schematically  a  damaged 
material  with  matrix  cracking  and  partial  debond¬ 
ing  near  the  crack  tips.  If  the  crack  spacing  (/) 
and  debond  length  (d)  are  given,  the  transverse 
shear  modulus  G23  can  be  determined  for  the 
given  applied  stress  by  using  a  micromechanics 
approach  based  on  Eshelby’s  equivalent  inclusion 
method  (see  Eshelby,  1957;  Mori  and  Tanaka, 
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1973;  Takahashi  and  Chou,  1988).  An  isotropic 
matrix  with  transversely  isotropic  ellipsoidal  in¬ 
clusions  is  considered. 

Assuming  small  displacements,  the  relation¬ 
ship  between  stress  (<7}y)  and  strain  (e,y)  tensors 
for  a  generally  anisotropic  material  is  expressed 
in  the  form 

~  ^ijkl^kl’  (^) 

or 

^  ijkl^kl' 

where  is  the  stiffness  tensor.  The  effective 
elastic  compliances  of  the  composite  5,^^/  are 
obtained  by  considering  the  work  done  under  the 
applied  stress  o;y: 

Sijki^^u^ki  =  S^kio-ijor^i  +  ,  (5) 

where  is  the  matrix  elastic  compliance  and 
e*  is  the  solution  of  the  following  simultaneous 
equations: 

+ (1  -  - «:,)] 

-  +  (1  - 

(6) 

where  and  are  the  stiffnesses  of  the 


matrix  and  fiber,  respectively,  and  is  known 

as  Eshelby’s  tensor. 

For  the  case  of  perfect  bonding  between  fiber 
and  matrix,  the  transverse  shear  modulus  G23  is 
obtained  by  substituting 

^23  ~  ^32  ~  1 ! 

all  other  a-,y  =  0,  into  Eq.  (5): 

0^3  GJt-Kf)  +  G23f(3-4..^  +  Kf) 

GJl  +  (3-4^,)Ff]  +  G,3K3-4..„)(l-K,)’  ^  ^ 

where  ^m>  ^23f  are  the  matrix  Poisson’s 
ratio,  matrix  shear  modulus,  and  transverse  shear 
modulus  of  the  fiber,  respectively.  Subscripts  2 
and  3  denote  the  transverse  in-plane  and  trans¬ 
verse  out-of-plane  directions,  respectively. 

For  the  case  of  complete  debonding,  the  trans¬ 
verse  shear  loading  is  replaced  by  the  equivalent 
tension  and  compression  at  45°  with  the  direction 
of  shear.  Using  the  same  equations  as  above  for 
this  case  we  obtain  shear  modulus  G23  as: 

G\^  (1-Kf)d> 


Fig.  2.  Composite  element  with  matrix  cracking  and  partial  debonding. 
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where 


where 


cl,  =  il+vj[\  +  2il-vJQ] 

X  (1  -  F,)-  +  [4R  +  (3  -  -  2v^v,,,)Q] 

X(l-K,)Kf  +  4(l-,.„) 

X  (1  —  >^i2f*^21f)(^ 

4(1  —  f^l2f’^21f)-^m 

(1  “  '^12f*^21f)'^m 


X  [(1  -  »'m^l2f)0  +  4(1  -  V„,I^i2f)] 
X(l-V,)V,R+^-^E,,VfQ. 


In  the  case  of  partial  debonding  along  the 
fiber  direction,  it  is  assumed  that  the  fiber-ma¬ 
trix  interface  around  the  circumference  of  the 
fiber  is  completely  debonded.  The  transverse 
shear  modulus  is  then  calculated  by  linear  inter¬ 
polation  of  the  moduli  for  the  cases  of  perfect 
bonding  and  complete  debonding,  i.e.. 


G23  =  G'’ 


I -2d  ,  2d 

23  ^  *■  ’ 


or 

Gn 


G« 


=  1  - 


1  - 


G’  ^ 


'23 


r'O 
^23  / 


a. 


(9) 


(10) 


This  approach  is  valid  when  the  crack  spacing 
is  sufficiently  smaller  than  the  characteristic  di¬ 
mension  of  the  ultrasonic  sensor  to  be  discussed 
below. 

The  ratio  of  debond  length  to  crack  spacing 
id/l)  is  obtained  from  a  modified  shear  lag  anal¬ 
ysis  model  proposed  by  Lee  and  Daniel  (1992). 


2d  1 

—  =  1 - log 

/  /3 


1+^ 

1-r 


(11) 


2F;. 


Et 


1)  £„o-j  +£,<r„ 


and 


V  = 


Ml 

^ - , 

'•f 


M=  y  H 


+ 


\E,  E^(v^-l)J’ 


H  = 


1  (3i7  +  l)('h  -  1) 

+ 


4Gf  3(77-1-  l)”^n 


(12) 


in  which  is  the  fiber  radius,  o-,  is  the  applied 
stress,  o^n,  is  the  axial  residual  stress  in  the  ma¬ 
trix,  £j5  is  the  interfacial  shear  strength,  and  £f, 
£„,,  Gf,  and  G„  are  the  Young’s  moduli  and 
shear  moduli  of  the  fiber  and  matrix,  respectively. 

Wooh  and  Daniel  (1991)  measured  the  five 
elastic  constants  of  a  transversely  isotropic  unidi¬ 
rectional  laminate  by  means  of  ultrasonic  wave 
propagation  tests.  Determination  of  Gnj  is  exper¬ 
imentally  much  simpler  than  that  of  the  other 
constants.  For  a  plane  wave,  the  displacements 
are  given  as 

«m=  exp[A:j(/t,x,-cr)],  (13) 

where  k  is  the  wave  number,  j  is  the  unit  imagi¬ 
nary  number,  n,  are  the  direction  cosines  of  the 
normal  to  the  wave  front,  is  the  displacement 
amplitude,  and  c  the  wavespeed.  Substituting  Eqs. 
(3)  and  (13)  into  the  equation  of  motion  yields  the 
eigenvalue  equation  and  noting  that  a  nontrivial 
solution  exists  when  the  determinant  vanishes, 
the  following  determinantal  equation  is  obtained 

\Ciikininj-  =  (14) 

where  5,^  is  the  Kronecker  delta  and  p  is  the 
density  of  the  material.  The  stress-strain  rela¬ 
tionship  for  orthotropic  materials  is  expressed  in 
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terms  of  nine  independent  elastic  constants,  that 
is: 


^11 

C,3 

0 

0 

0 

^22 

^12 

^22 

^23 

0 

0 

0 

^22 

T 

^23 

,  = 

^13 

0 

^23 

0 

^33 

0 

0 

C44 

0 

0 

0 

0 

^33 

723 

0 

0 

0 

0 

Qs 

0 

731 

0 

0 

0 

0 

0 

Ti2J 

(15) 

where  [C,y]  is  the  stiffnes  matrix,  subscripts  1,  2,  3 
denote  the  principal  material  axes,  and  are 
the  engineering  shear  strains.  Then,  Eq.  (14)  takes 
the  form 


where  the  so-called  Christoffel  stiffnesses  are 
given  by 


Ai 

+  215055, 

^22 

^  ^2^22 

+  /I3C44, 

As 

~  ^2^44 

+  215033, 

(17) 

As 

=  «1«2(<^12  +  Q6)» 

^23 

=  n2fi  2(023  +  C44), 

r,3 

=  «3«i(A3  +  Q5)- 

■ 

By  setting  the  directional  cosines  for  a  wave 
motion  traveling  in  the  thickness  direction  Xy, 
i.e., 

ni  =  /i,  =  0  and  n^=\,  (18) 


^11  ~ 


det 


^2 

^3 


rn 

pc~ 

As 


As  p^~ 


(16) 


three  distinct  solutions  for  c  are  obtained.  Of 
these,  the  solution  of  interest  is  the  one  related 
to  the  transverse  shear  modulus,  that  is 

G23  =  C,,=pcl„  (19) 


A 


Grip 


Vi 


Transmitter 


Strain 


iJ 


Receiver 


Load  controller 


Signal  conditioner 


Digital  oscilloscope 


Arbitrary  function  generator 


Power  amplifier 


Fig.  3.  Schematic  block  diagram  of  the  ultrasonic  system. 
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where  is  the  velocity  of  the  shear  wave  polar¬ 
ized  in  the  plane  transverse  to  the  fiber  direction. 
Fiber-matrbc  debonding  with  a  surface  normal 
transverse  to  the  wave  propagation  direction  af¬ 
fects  the  transverse  shear  modulus  degradation 
which  is  related  to  the  corresponding  wavespeed. 


3.  Experimental  procedure 

The  material  investigated  was  a  unidirectional 
SiC/CAS,  calcium-aluminosilicate  glass  ceramic 
reinforced  with  silicon-carbide  fibers,  manufac¬ 
tured  by  Coming  Glass  Works.  The  fiber  is  sili¬ 
con-carbide  yarn  known  as  Nicalon  (Nippon  Car¬ 
bon  Co.)  of  15  fim  average  diameter.  The  com¬ 
posite  was  obtained  in  the  form  of  16-ply,  3.048 
mm  (0.12  in)  thick  unidirectional  plates. 

Unidirectional  coupons,  3.05  mm  (0.12  in) 
thick,  17.27  mm  (0.68  in)  wide,  and  15.24  cm  (6 
in)  long,  were  prepared.  They  were  tabbed  with 
2.54  cm  (1  in)  long  glass/ epoxy  tabs  and  instru¬ 
mented  with  strain  gages  on  both  sides.  The 
specimens  were  loaded  in  a  servohydraulic  testing 
machine  while  applied'  stress  and  strains  (both 
longitudinal  and  transverse  strains)  were  moni¬ 
tored  and  recorded.  Meanwhile,  a  pair  of  ultra¬ 
sonic  shear  wave  transducers  were  tightly  at¬ 
tached  to  the  surfaces  of  the  sample  for  ultra¬ 
sonic  characterization.  These  transducers  were 


set  up  in  through-transmission  mode  and  they 
were  aligned  to  generate  a  shear  wave  polarized 
perpendicular  to  the  fiber  direction.  A  micro¬ 
computer  was  used  to  control  the  entire  process 
of  load  control  and  data  acquisition.  A  10  MHz 
toneburst  signal  generated  by  an  arbitrary  func¬ 
tion  generator  (LeCroy  9100)  excited  the  trans¬ 
mitter  after  passing  through  an  RF  power  ampli¬ 
fier  (ENI  325LA).  Then,  the  transmitted  multiple 
echoes  were  recorded  by  a  digital  oscilloscope 
(Tektronix  RTD710)  mnning  at  a  sampling  rate 
of  200  MHz  (5  ns  intervals).  Figure  3  shows  a 
schematic  diagram  of  the  ultrasonic  test  system 
used.  As  soon  as  all  mechanical  and  ultrasonic 
measurements  were  completed,  the  specimen  was 
immediately  loaded  to  the  next  step  and  this 
procedure  was  repeated  until  global  failure  of  the 
specimen. 


4.  Wavespeed  measurements 

Wavespeeds  were  measured  from  the  wave 
traces  of  the  transmitted  and  multiple  echoes  in 
the  material  as  illustrated  in  Fig.  4.  Figure  5  is  a 
typical  as-obtained  waveform.  The  wavespeed  is 
simply 

c  =  2h/At,  (20) 


Transmitter  Specimen  Receiver 


Gates:  1  2  3 

rn  r~i  rn 


Before  loading 


Fig.  4.  Transmitted  signals  and  multiple  echoes  for  measurement  of  wavespeed  and  attenuation. 
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Fig.  5.  Illustrative  example  of  cepstrum  analysis  for  measuring 
wavespeed:  (a)  as-obtained  waveform;  (b)  magnitude  spec¬ 
trum;  and  (c)  magnitude  cepstrum. 


where  Af  is  the  time  interval  between  echoes  and 
h  is  the  specimen  thickness.  The  time  interval  can 
be  directly  obtained  in  the  time  domain  by  mea¬ 
suring  the  time  interval  between  two  correspond¬ 
ing  peaks  of  two  consecutive  echoes  because  the 
signal  used  is  a  single  frequency  toneburst.  How¬ 
ever,  it  is  not  desirable  to  measure  the  wavespeeds 
in  this  manner  since  it  is  sometimes  difficult  to 
locate  the  right  peaks  and  the  complex  algorithm 


for  automated  testing  increases  the  chance  of 
introducing  errors. 

An  alternative  approach  is  to  determine  the 
wavespeed  based  on  ultrasonic  interference  spec¬ 
troscopy  (see  Chang  et  al.,  1974;  Wooh  and 
Daniel,  1991).  In  this  method,  the  obtained  wave¬ 
forms  fit)  are  expressed  in  the  frequency  do¬ 
mains  as 

F{a))  =  (  fir)  exp(-ymT)  dr,  '(€1) 

— QO 

by  means  of  fast  Fourier  transform  (FFT)  (see 
Brigham,  1988).  The  relatively  low  periodicity  of 
the  echoes  compared  to  the  frequency  of  the 
ultrasonic  wave  produces  well  defined  anti-reso¬ 
nance  dips  and  peaks  in  the  frequency  spectrum 
as  shown  in  Fig.  5b.  The  period  of  these  dips  or 
peaks  directly  gives  the  time  interval  Ar.  One  of 
the  advantages  of  this  method  is  that  the  dips 
appear  even  when  the  echoes  interfere  with  each 
other,  that  is,  this  technique  is  useful  even  for 
very  thin  laminates.  However,  because  of  the 
many  different  situations,  the  algorithm  of  search 
for  the  local  minima  is  computationally  complex. 
The  method  of  fails  in  a  noisy  environment,  espe¬ 
cially  in  a  case  when  the  specimen  has  severe 
internal  damage  and  the  signal-to-noise  ratio  be¬ 
comes  low. 

To  solve  this  problem,  the  so-called  cepstrum 
analysis  was  used  (see  Childers  et  al.,  1977;  Wooh 
and  Daniel,  1992).  The  above-mentioned  spec¬ 
trum  is  considered  as  another  time  series.  Taking 
the  FFT  of  the  magnitude  spectrum  yields  a 
spectrum  of  the  spectrum  called  magnitude  cep¬ 
strum  which  is  back  in  the  time  domain: 

g(t)  =  f  |F(ft))|  exp(-;aj/)  dr.  (22) 

•'  —00 

The  time  in  the  cepstrum  is  called  quefrency  and 
since  the  interferences  in  the  magnitude  spec¬ 
trum  are  periodic,  the  maximum  amplitude  oc¬ 
curs  at  the  quefrency  f  =  A/  in  the  magnitude 
spectrum  as  illustrated  in  Fig.  5c.  The  algorithm 
for  determining  a  time  delay  becomes  much  sim¬ 
pler  because  one  needs  to  simply  find  a  maximum 
amplitude  and  its  absolute  location  in  the  que- 
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frency  domain.  Wavespeeds  are  then  calculated 
from  Eq.  (20). 

A  drawback  of  this  method  is  that  it  is  compu¬ 
tation  intensive  because  one  needs  to  perform 
the  Fourier  transformation  twice  for  each  loading 
step.  Fortunately,  modem  computers  are  capable 
of  handling  such  a  task  in  real  time. 


5.  Attenuation  measurements 

Attenuation  is  another  important  parameter 
indicating  damage  in  the  material.  Although  we 
do  not  have  an  analytical  solution  for  attenuation, 
it  can  be  measured  and  correlated  with  the  inter¬ 
nal  damage.  In  contrast  to  the  transverse 
wavespeed,  the  attenuation  is  sensitive  to  matrix 
cracking  as  well.  The  amplitudes  of  the  first  three 
echoes  in  Fig.  4  are: 

=  kAqT^qTqq{}, 

"^2  ~ (23) 

A^  =  '^^O^AB^BC-^BC^BA'^^’ 

where  k  is  the  sensitivity  of  the  receiving  trans¬ 
ducer,  Aq  the  amplitude' of  the  incident  signal, 
T^b,  Tbc,  ^ba>  ^bc  ^re  the  transmission  and 
reflection  coefficients  at  the  front  and  back  sur¬ 
faces  of  the  specimen,  respectively.  The  signal 
loss  is  then  expressed  by  solving  the  above 
equations: 


The  attenuation  can  be  expressed  as  the  sum  of 
the  initial  attenuation  (ofp;  no  damage)  and  the 
attenuation  due  to  damage  (a^): 


a  =  «o  +  «d-  (25) 

Therefore, 

20  [  n\ 

a,  =  a-ao=--log^-J 
20 

=  -  y  logio^,  (26) 


where  the  relative  signal  loss  q  is 

A,  I  ^,/Ai 

^  (^i)o  {(^2/A)oI 

(  ^3/^2  I  A, /A, 

\(^3/^2)oj  \(^3/^l)o/  ^ 

Subscript  zero  means  initial  reference  stage  be¬ 
fore  loading.  One  can  select  any  two  echoes  to 
calculate  the  attenuation,  but  it  is  more  desirable 
to  pick  the  echoes  that  correspond  to  the  longest 
travel  in  order  to  minimize  measurement  errors. 

Amplitudes  of  the  individual  wavelets  were 
measured  by  introducing  an  imaginary  gate  on 
the  waveforms  and  finding  peak-to-peak  values 
within  the  gate  (Fig.  4).  The  location  of  the  gate 
is  recalculated  by  the  wavespeed  measurement  so 
that  the  gate  follows  the  corresponding  wavelets. 


6.  Results  and  discussion 

Some  properties  of  the  matrix  and  fiber  con¬ 
stituents  obtained  from  the  literature  and  from 
tests  are  tabulated  in  Table  1  (Larsen  and  Adams, 
1988).  The  unidirectional  composite  was  mechan¬ 
ically  characterized  and  the  results  are  shown  in 
Table  2. 

As  discussed  before,  the  method  used  here  is 
valid  when  the  debond  length  is  sufficiently 
smaller  than  the  characteristic  dimension  of  the 
ultrasonic  transducer.  In  the  present  case,  the 
diameter  of  the  transducer  used  was  25.4  mm  (1 


Table  1 

Constituent  material  properties 

Property  CAS  SiC 

_ _  matrix  fiber 

Maximum  use  temperature,  1350  (2460)  1300  (2370) 

°C  (°F) 

Fiber  diameter,  p.m  -  15 

Density,  g/cm^  2.8  2.6 

Coefficient  of  thermal  expansion,  5.0  (2.8)  3.2  (1.8) 

10“V‘'C  (10"V°F) 

Elastic  modulus,  GPa  (Msi)  98  (14.2)  170  (280) 

Tensile  strength,  MPa  (ksi)  124  (18)  1930  (280) 

(flexural) 
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Table  2 

Measured  properties  of  SiC/CAS  unidirectional  composite 


Property 

Value 

Fiber  volume  ratio, 

0.39 

Ply  thickness,  mm  (in) 

0.38  (0.015) 

Longitudinal  modulus,  GPa  (Msi) 

125(18.1) 

Transverse  modulus,  £'2»  (Msi) 

112(16.2) 

In-plane  shear  modulus,  MPa  (Msi) 

52  (7.5) 

Major  Poisson’s  ratio, 

0.18 

Longitudinal  tensile  strength,  MPa  (ksi) 

393  (57) 

Transverse  tensile  strength,  F2t,  MPa  (ksi) 

55  (8) 

Longitudinal  ultimate  tensile  strain, 

0.0084 

Transverse  ultimate  tensile  strain, 

0.0005 

in)  which  is  approximately  160  times  the  mea- 
sured  minimum  crack  spacing  (160  jim).  Figure  6 
shows  the  predicted  and  measured  variation  of 
the  transverse  shear  modulus  G^i  due  to  induced 
debonding  during  loading.  It  shows  qualitatively 
that  the  modulus  is  degraded  as  debonding  grows. 
However,  from  the  quantitative  point  view,  the 
experimental  curve  shows  some  deviation  from 
the  prediction.  The  experimentally  measured  rate 
of  degradation  is  more  uniform  than  the  pre¬ 
dicted  one.  This  disagreement  could  be  attributed 
to  the  neglect  of  the  frictional  stress  between  the 
fiber  and  matrix.  The  friction  prevents  the  sud¬ 
den  propagation  of  debonding  immediately  after 
initiation  resulting  in  the  less  steep  curve  in  this 
region.  The  agreement  is  better  in  the  last  part  of 


Fig.  6.  Change  in  transverse  shear  modulus  due  to  debonding 
during  loading. 


Fig.  7.  Change  in  normalized  signal  loss  due  to  induced 
damage  during  loading  with  superimposed  stress-strain  curve. 

the  curve  corresponding  to  more  extensive 
debonding. 

Figure  7  shows  the  change  in  normalized  ultra¬ 
sonic  signal  loss  due  to  induced  damage  during 
loading  as  a  function  of  applied  strain.  The  atten¬ 
uation  shows  a  good  qualitative  correlation  with 
previously  observed  failure  mechanisms  and  with 
the  associated  stress-strain  curve.  The  point  of 
initial  change  in  signal  loss  coincides  exactly  with 
the  proportional  limit. 

In  summary,  ultrasonic  methods  were  devel¬ 
oped  for  real-time  monitoring  of  damage  evolu¬ 
tion  in  ceramic  composites  under  tensile  loading. 
The  transverse  shear  modulus  calculated  from 
wavespeed  measurements  can  be  used  as  a  sensi¬ 
tive  indicator  of  fiber- matrix  debonding.  Ultra¬ 
sonic  attenuation  shows  good  qualitative  agree¬ 
ment  with  damage  and  the  combination  of  atten¬ 
uation  and  wavespeed  would  be  a  useful  method 
to  study  failure  mechanisms.  Future  studies 
should  include  the  effects  of  residual  stresses  and 
the  friction  between  fiber  and  matrix  for  a  better 
prediction  of  material  behavior. 
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Abstract— A  reduced  interphase  stiffness  of  a  single  fiber  (the  dissimilar  fiber)  in  a  unidirectionally 
reinforced  composite  gives  rise  to  stress  deviations  in  its  own  interphase,  as  well  as  in  the  interphases 
of  neighboring  fibers,  relative  to  the  stresses  in  a  perfect  composite.  For  transverse  loading  and  an 
arbitrary  cross-sectional  distribution  of  the  fibers,  a  general  method  is  presented  to  calculate  the  e 
stress  deviations,  based  on  solutions  by  the  boundary  element  method  '"'®8ra 

for  the  dissimilar  fiber  and  neighboring  fibers.  In  this  method  nearest 

taken  into  account.  The  interphases  are  represented  by  the  spring  layer  model.  Det.nled  ^ 

leLts  are  presented  for  the  special  case  of  a  hexagonal  array  composite.  Results  are  compared  for 
ralculafionrtaking  into  account  nearest  and  next-nearest  neighbors.  Stresses  at  the  inatnx  sides  of 
the  interphases  and  energy  densities  in  the  interphases  have  been  calculated  for  the  dissimilar  fiber 
Ind  ?or  S  nextneares ^neighbors.  These  stresses  have  also  been  obtained  tor  the  case  that  the 
dissimilar  liber  has  interphase  Daws. 


INTRODUCTION 

The  overall  mecfianical  properties  and  the  strength  of  fiber-reinforced  composites  are 
significantly  affected  by  the  properties  of  very  thin  layers  at  the  fiber-matnx  interfaces  i.e. 
by  fiber-matrix  interphases.  Several  authors  have  investigated  the  effects  of  interphase 
compliance  and  interphase  flaws  on  the  effective  elastic  constants  of  the  composite  material. 
In  these  studies  interphases  have  generally  been  represented  by  the  spring-layer  modeh  In 
this  model  the  interphases  are  treated  as  very  thin  zones  of  unspecified  thickness.  The  radial 
and  tangential  tractions  are  continuous  across  the  interphase,  but  the  displacements  inay 
be  discontinuous  between  fiber  and  matrix  due  to  the  presence  of  the  mterphase.  It  has 
generally  been  assumed  that  the  tractions  are  proportional  to  the  corresponding  dis¬ 
placement  discontinuities.  The  proportionality  constants  characterize  the  stiffness  of  the 
interphase.  This  “spring-layer  model”  has  been  employed  by  Benventste 

(1985),  Steifand  Hoysan  (1987),  Hashin  (1990)  and  Achenbach  and  Zhu  (1989).  Benveniste 
(1985)’ calculated  the  effective  modulus  of  a  composite  reinforced  by  spherical  particles 
Uich  are  not  well  bonded  to  the  matrix.  Steifand  Hoysan  (1987)  used  an  energy  method 
for  calculating  the  longitudinal  stiffness  of  aligned  short-fiber  composites  with  imperfectly 
bonded  interfaces.  Hashin  (1990,  1991)  used  the  spring  interphase  model  in  hts  analysis  of 
the  thermoelastic  behavior  of  a  fiber  reinforced  composite  as  well  as  that  of  a  particula  e 
composite  He  (1992)  applied  extremum  principles  of  the  theory  of  elasticity  to  composite 
bodies  to  obtain  simple  bounds  for  the  effective  elastic  properties  of  tw^phase  materials 
with  imperfect  interfaces.  Jasiuk  et  al.  (1992)  have  investigated  the  eflect  of  a  sliding 
interface  on  the  elastic  properties  of  composites  with  randomly  distributed  circular  and 
spherical  rigid  inclusions.  Recently,  Gosz  et  al.  (1992)  obtained  the  transverse  constitutive 
response  of  a  hexagonal  array  composite  using  a  combined  analytical  and  numerical 

'"^^'^The  effect  of  interphase  flaws  and  radial  matrix  cracks  on  the  overall  mechamcal 
properties  has  been  studied  by  Achenbach  and  Zhu  (1989,  1990)  and  Achenbach  and  Choi 
(1991)  These  investigations  have  been  carried  out  for  periodic  arrays  of  the  fibers,  such  a 
rectangular  and  hexagonal  arrays,  and  it  has  been  assumed  that  interface  flaws  and  matrix 
cracks  have  the  same  periodic  distribution  as  the  fibers.  By  virtue  of  these  geometrica 
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simplifications  it  was  possible  to  consider  a  basic  cell  for  detailed  calculations  of  the  fields 
of  stress  and  deformation. 

In  this  paper  the  more  general  case  is  considered  of  a  single  fiber  which  has  either  a 
smaller  interphase  stiffness,  or  whose  interphase  contains  a  defect.  All  the  other  fibers  have 
the  same  interphase  properties.  The  stress  deviations  in  the  interphases  generated  by  the 
presence  of  the  dissimilar  fiber  are  calculated  and  relevant  stresses  and  deformation  energy 
densities  in  the  interphases  of  the  composite  with  the  dissimilar  fiber  are  compared  with  the 
corresponding  quantities  in  the  perfect  composite.  This  work  generalizes  an  earlier  study 
by  Zhu  and  Achenbach  (1991).  The  approach  used  in  this  paper  is  applicable  for  linearly 
elastic  behavior  of  the  fibers  and  the  matrix. 


FORMULATION 

Figure  1  shows  a  cross-sectional  view  of  a  unidirectionally  reinforced  fiber-composite 
in  which  all  fibers  are  of  the  same  radius  a  and  have  the  same  interphase  properties,  except 
one:  the  dissimilar  fiber.  This  fiber  has  different  interphase  stiffness  properties  including 
the  possibility  that  all  or  part  of  the  interphase  may  have  zero  stiffness,  i.e.  the  fiber  may 
be  completely  or  partially  debonded.  It  is  assumed  that  at  a  remote  location  the  composite 
is  subjected  to  uniform  normal  stresses  Pa^  and  Qgq  in  the  x  and  y  directions,  respectively. 
The  fibers  are  labelled  by  the  index  ‘Vy”,  where  the  index  “0”  is  used  to  identify  the  dissimilar 
fiber,  and  where  roughly  speaking  the  larger  the  value  of  the  larger  the  distance  from  the 
dissimilar  fiber.  The  circumference  of  fiber  q  is  denoted  by  T^.  A  large  contour  Too  is  the 
boundary  of  an  area  of  interest  of  the  cross  section  shown  in  Fig.  1.  Inside  F^o  there  are 
Af-f  1  fibers.  Thus  ^  =  0,  I,  2, . . . ,  A. 

Following  earlier  papers  by  Achenbach  and  Zhu  (1989,  1990),  a  very  thin  compliant 
interphase  between  fibers  and  matrix  is  modeled  by  a  distribution  of  mechanical  springs 
across  a  representative  interface.  This  model  implies  that  the  tractions  are  continuous,  but 
the  displacements  may  be  discontinuous  across  this  interface.  For  linear  behavior  in  the 
interphase,  the  relations  between  the  relevant  traction  and  displacement  components  may 
then  be  written  as  (at  r  =  a) 


(a)  (b) 

Fig.  1.  (a)  Cross-sectional  view  of  the  composite  with  a  fiber  with  dissimilar  interphase  stiffness, 

(b)  configuration  of  a  fiber. 
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where  t,  and  te  are  the  interface  tractions  in  the  radial  and  circumferential  directions, 
respectively.  Here  and  in  the  sequel,  quantities  with  upper  index  “m”  and  “f”  are  defined 
in  the  matrix  and  the  fiber  regions,  respectively.  The  addition  of  eqn  (lb)  assures  that  the 
model  will  not  allow  a  physically  unrealistic  radial  overlap  of  the  matrix  and  fiber  materials 
across  the  interface.  The  constants  and  kg  are  proportionality  constants  which  define  the 
mechanical  properties  of  the  interphase. 

The  conditions  (1)  include  the  case  of  perfect  bonding  (/:,  =  ko=  oo)  when  the  tractions 
and  displacements  are  continuous  and  the  case  of  total  debonding  (k^  =  ko  =  0)  when  the 
tractions  vanish.  In  the  latter  case,  it  should  be  noted  that  the  pure  sliding  condition 
(k,  =  oo,ki  =  0)  is  assumed  to  exist  on  the  interface  when  the  radial  stresses  are  compressive. 
If  there  exist  interphase  defects  like  open  cracks,  say  over  a  region  -0,<0<6,<  n/2  and 
-0^  +  n  <6  <  0^+n,  then  in  the  domain  0  <  0  ^  7c/2  the  conditions  (1)  are  valid  for  r  =  a, 
0e  <  0  <  7t/2,  while  for  r  =  n  and  0  <  0  <  0,;,  we  should  have  interphase  crack  conditions 
defined  by 


/,"*  =  =  /S'  =  tS  =  0. 


(2) 


It  is  also  noted  that  for  an  interphase  crack  the  ligament  at  the  tip  of  the  disbond 
undergoes  a  finite  stretch  when  in  tension,  and  consequently  the  tractions  remain  bounded. 
Hence  the  usual  problems  of  violently  oscillating  singularities  (see  Williams,  1959)  that  are 
associated  with  crack-tip  fields  for  a  crack  in  a  perfectly  bonded  interface,  do  not  occur 
for  the  spring-layer  model.  This  conclusion  follows  immediately  from  eqn  (1)  and  the 
boundedness  of  the  displacements. 

For  the  generation  of  interphase  cracks,  as  well  as  their  propagation  and  arrest,  it  is 
feasible  to  use  a  critical  stress,  critical  strain,  or  critical  strain  energy  density  criterion, 
because  in  the  spring-layer  model  these  quantities  are  well  defined  near  the  tip  of  an 
interphase  crack.  In  this  paper  we  will  employ  an  energy  density  criterion,  since  it  combines 
information  on  the  tensile  and  shear  stresses  in  the  interphase.  For  the  interphase  model 
defined  by  eqns  (I),  the  strain  energy  per  unit  interphase  area  is  easy  to  calculate.  We  have 


2K  2ko 


(3) 


It  should  be  noted  here  that  ll  is  included  in  U  only  when  t[  is  positive  (tension).  It  is 
assumed  that  compressive  values  of  t\  do  not  affect  the  integrity  of  the  interphase. 

By  substituting  k  into  eqn  (3)  the  normalized  form  oi  V,V  =  Ufi’^Kala),  becomes 


(J  ^  OeKy 

2k,a/n"'  2koaUi'” 

where  fi”'  is  the  shear  modulus  of  the  matrix  material,  and  tr#  is  the  applied  stress.  It  is 
reasonable  to  assume  that  the  interphase  will  break  and  form  an  interphase  crack  in  a 
region  of  positive  radial  stress  when 


U  >  f/“.  (5) 

When  U  =  f/",  the  applied  stress  Oo  reaches  a  critical  value  tr"  which  can  be  related  to  U" 
by  eqn  (4)  as 


a  n\ 


(6) 


An  analysis  of  the  initiation  and  propagation  of  matrix  cracks  must  also  be  based  on 
an  appropriate  criterion.  For  a  perfect  composite  subjected  to  tensile  stresses,  numerical 
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results  show,  in  agreement  with  physical  intuition,  see  e.g.  Achenbach  and  Zhu  (1989, 
1990),  that  the  circumferential  tensile  stress  at  the  fiber-matrix  interphase  is  the  largest 
tensile  stress  component  in  the  matrix  material.  As  a  crack  initiation  criterion  we  choose 

oo  >  (7) 


where  <t"  is  a  critical  stress  of  the  matrix  material  for  matrix  crack  initiation. 

The  displacement  and  traction  components  for  the  case  that  no  dissimilar  fiber  is 
present  (perfect  composite)  are  denoted  by  and  /f'‘(x)  for  the  fibers,  and  and 

t,“'’(jc)  for  the  matrix  material.  The  deviations  of  these  quantities  due  to  the  presence  of  a 
dissimilar  fiber  are  denoted  by  flf(x),  ff(x),  ft^ix)  andyVW-  The  total  displacement  and 
traction  fields  in  the  presence  of  a  dissimilar  fiber  may  then  be  written  as 


i/f(.v)  =  u[>’(.x)  +  i7f(A:) 

(8a) 

m;"(jc)  =  m;"‘’(x:)  +  «7"(x) 

(8b) 

t%x)  =  /?’(A:)  +  r[(.x) 

(8c) 

tr(A:)  = /r(^)  +  tT(A:). 

(8d) 

The  interface  conditions  corresponding  to  the  case  of  a  different  interface  constant  are 
obtained  by  replacing  in  eqns  (la),  (lb)  and  (Ic)  the  total  tractions  t”,  t^,  t[,  and  by  the 
superpositions  (8c)  and  (8d).  On  the  interphase  To)  we  then  have; 

-if  = 

(9a) 

-fir  =  koiu'o-rfu) 

(9b) 

where  (r,  0)  is  a  local  polar  coordinate  located  at  the  center  of  each  fiber  and 
the  interfacial  boundary  of  the  fiber  q.  On  the  interphase  To  we  have : 

represents 

-if  =  -(k-kMir-iif)+kX{ef-ii) 

(9c) 

-ie  =  -{ko-h)iiir-tih+koiito-u[). 

(9d) 

In  eqns  (9c)  and  (9d),  kr  and  ks  are  the  interphase  constants  for  the  dissimilar  fiber. 

If  the  interphase  Vq  contains  an  interphase  crack  in  the  range  of  —d^<6<0c<  njl 
and  n  —  0c<  9  <  then  eqns  (9c  and  9d)  become 

-if  =  -kr(lC-llf)  =  tfP 

(9e) 

a‘ 

II 

T 

1 

II 

1 

(9f) 

For  all  interphases  F,,  including  Fq,  the  continuity  of  the  tractions  across  the  interphase  F, 
are  satisfied  by 

—  f =  f 

If  » r 

(10a) 

_  7m  __  7f 
—  iQ  —  If). 

(10b) 

By  appropriate  coordinate  transformations  the  interphase  conditions  (la),  (lb)  and 
(Ic)  can  be  expressed  in  terms  of  Cartesian  components  of  tractions  and  displacements 
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which  will  be  used  in  the  actual  computations.  The  interphase  conditions  on  the  interphase 
Tq  (#  To)  given  by  eqns  (9a,  9b)  may  be  rewritten  as 


nr. 


(11a) 


f" 

—  aC2  cos  —  +  aC,  sin  0 


jm 

‘y 


sin  0{u'^  —  fij.) +COS  flOTjr  —  =  0 


(11b) 


where  has  been  introduced  to  enhance  the  accuracy  of  the  solutions  of  the  simul¬ 

taneous  equations  in  the  numerical  calculations  using  the  BEM.  The  dimensionless  com¬ 
pliant  constants  are  defined  as 


1  /("’  1  u”' 

C,  —  =  —  and  C2  =  —  =  — . 
k\  cikr  ki  akf) 


(12a,b) 


Similarly  the  interphase  conditions  (9c,  9d)  on  the  boundary  of  the  dissimilar  fiber. 
To,  can  be  expressed  as 


f 

aC|  cos0-^ -t-aC,  sin0—  +cos  0(/7'^  — H()  +  sin  0(0^— Oy) 

/j”  ft"' 


^mp 

“}-a(C|  “Cj)  [  cosO— — f-sin  9  )  =  0  (13a) 


—aC2  sin  0—  ^aCj  cos  0  —  —  sinfl(r/™  — »v)+cos  — /”/[,) 

/i'”  /i"’ 


-a(C2-C2)  ( sin0“  -COS0  — )  =  0  (13b) 


where 


1  //■"  -  1  /i”' 

£|  aH,  ~  Ic2~  aica ' 


(14a,b) 


If  the  interphase  of  To  contains  an  interphase  crack  in  the  range  o(  —6y<6  <dy<  njl 
and  n  —  6^<  0  <n  +  6^  then  eqns  (13a,  13b)  become 


(I5a,b) 


If  the  dissimilar  interphase  To  has  compressive  radial  tractions  over  one  or  more 
elements  on  To  then  to  prevent  the  interphase  from  overlapping,  the  eqns  (13a,  13b)  become 


Wx  =  u*  and  mJ"  =  i7j,  on  To  if <  0. 
Continuity  of  tractions  across  the  interphase  F,  including  Fo  is  satisfied  as 


(16a,b) 


and  ty  = -f,. 


(17a,b) 
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BIEs  for  unidirectional  composite  with  one  dissimilar  fiber 

In  earlier  papers,  see  e.g,  Achenbach  and  Zhu  (1989,  1990),  a  boundary  integral 
equation  has  been  derived  for  the  displacements  and  tractions  on  the  circumference  of  a 
fiber  inside  a  basic  cell  of  the  composite.  Here  we  write  the  analogous  boundary  integral 
equation  for  fiber  q. 

f/;"(A:,^)/7(jc)dr(jr)-f  7-(x,^)i/f(jr)dr(x),  xeT^^er,  (18a) 

r+r^  Jr+r^ 


where 

N 

I  r, 
,  =  0 

and 


1 _  ■ 

1-v'")  . 


(3-4v'")ln 


(18b) 


(19a) 


U^(x,OSj,+fi'^ 


dxj 


Ut(x,0  n*.(jc). 


(19b) 


A  similar  blE  for  the  solution  Hr'’(iJ)  in  the  matrix  of  the  perfect  composite  can  be  written 
as 


«r(?)  = 


c/s’(^.«)'r(-'^)<injc) 


r+r„ 


-f 

Jr+r,, 


7’y(jr.«)jr(jc)dr(x),  jcer.^er,.  (20) 


Subtracting  eqn  (20)  from  eqn  (18)  gives  a  BIE  for  ilTiO  ns 


\aT(lt)  = 


u^(x,OtT(x)mx)- 


T^(x,  |)Hf(jf)  dr(A:) 


+ 


f/^(;r.?)r7(x)dr(ji:) 


-[  T^(x,^X 

Jr„ 


(jc)dr(jr),  jrer.fer,.  (21) 


Since  the  tractions  corresponding  to  rC(x)  are  self-equilibrating  with  respect  to 
the  center  of  fiber  zero,  it  follows  from  Saint-Venant’s  principle  that  i7f(jc),  m™(ji^),  f/(jf)  nnd 
f["  (jc)  decrease  as  the  field  point  x  moves  away  from  the  center  of  the  dissimilar  fiber. 
Hence,  the  contour  can  be  appropriately  chosen  so  far  from  the  origin  that  the  integrals 
over  the  contour  r^,  in  eqn  (21)  come  to  vanish.  As  an  additional  approximation  it  may  be 
assumed  that  only  a  small  number  of  neighboring  fibers  need  be  included  in  F. 

The  corresponding  BIE  for  the  fiber  q  may  be  written  as 


C/f(jr,|)/j(jr)dr(.t)- 


T'iiix,^Wj{x)AT(x), 


(22) 


Similarly,  substituting  eqn  (8a)  into  eqn  (22)  yields 
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f/f(jf,?)/j(jc)dr(x)-  rf(x,^)i7}(x)dr(j:) 

Jr,  Jr, 

/* 

Ul(x,^)l'iHx)dr(x)-  rf(jr,^)»^-{;r)dr(jt),  AC.^er,.  (23) 

Jr,  Jr, 

But  we  know  that  the  terms  in  the  second  line  of  eqn  (23)  vanish  for  all  fibers  in  the 
composite  without  a  dissimilar  fiber,  so  we  have 

1  /* 

f/f>(x,?)f;(.t)dr(.t)-  n(ji-,?)i7Kx)dr(jr),  x,^6r,.  (24) 

For  plane  strain  the  circumferential  stress  crjf  along  the  matrix  side  of  the  interphase 
may  be  calculated  by 


where  du^jdO  will  be  obtained  by  a  numerical  differentiation  method. 

Example :  hexagonal-array  composite 

Equations  (l)-(25)  are  valid  for  any  distribution  of  the  fibers  over  the  cross  section. 
We  will  now,  however,  consider  the  special  case  of  a  hexagonal-array  composite  as  shown 
in  Fig.  2.  For  the  case  without  a  dissimilar  fiber  the  interphase  fields  have  been  analysed  by 
Achenbach  and  Zhu  (1989,  1990). 


Fig.  2.  Configuration  for  the  hexagonal  array :  (a)  the  dissimilar  fiber  0  has  lower  interphase  stiffness. 
All  other  fibers  have  the  same  interphase  stiflfnesses,  (b)  shows  the  domain  for  numerical  calculations, 
and  (c)  the  dissimilar  fiber  0  has  interphase  debonding. 
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The  centers  of  the  regular  fibers  are  located  at  the  following  points 


(^,v) 


|/|,|/|  =  0, 1,2,3, ...  and 


/-hy  =  even  integer 


(26) 


where  b  is  the  distance  between  the  centers  of  two  adjacent  fibers,  and  |/|  +  [/|  0.  A 

Cartesian  coordinate  system  (  v,  y)  has  its  origin  at  the  center  of  the  dissimilar  fiber. 

Because  of  symmetry  with  respect  to  the  jc-  and  ;'-axes  we  only  need  to  consider  the 
fibers  in  the  first  quadrant,  as  shown  in  Fig.  2.  The  nearest  neighbors  to  the  dissimilar  fiber 
are  denoted  fibers  1  and  2,  and  the  next  nearest  fibers  are  3,  4,  5  and  6. 

The  symmetry  conditions  on  the  sides  of  the  first  quadrant  make  it  possible  to  limit 
the  boundary  element  calculations  to  that  quadrant  only.  The  fundamental  solutions  which 
satisfy  the  symmetry  conditions  =  /^,  =  0  at  jc  =  0  and  iiy  =4  =  0  at  y  =  0  can  be 
constructed  from  the  full-space  fundamental  solutions  given  by  eqns  (19a)-(19b).  For  a 
load  applied  at  (a%  jO  the  construction  is  easily  achieved  by  placing  corresponding  loads  at 
the  image  points  of  (a,  v)  in  the  other  three  quadrants,  as  shown  in  Fig,  3.  In  terms  of 
Uij(x,^)  and  ^^(a,  given  by  eqns  (19a)-(19b)  (superscript  “m”  can  be  disregarded),  the 
first  quadrant  fundamental  solutions  C9(a,  ^)and  T9(a,  ^)  may  be  written  as 


U^(x,lt)  =  U,j(x,l^)-\-aUijix,  +  -^)+yUij(x,l),  x,^e  1st  Quadrant  (27a) 

T^(x,4)  =  T;j(x,!^)  +  aT,j(x,  +  x,^e  1st  Quadrant  (27b) 

where  the  coefficients  a,  ft  and  y  take  values,  oi  =  ft=—y=—\  for  /=1  and 
_a  “  ^  y  =  —  1  for  /  =  2  and  |  denotes  the  conjugate  point  of 

Morp'detailed  expressions  for  f/9(A,  ^)  and  T9(a,  ^)  are  given  in  the  Appendix. 

Boundary  element  method 

Equations  (18)~-(24)  have  been  solved  numerically  by  an  application  of  the  boundary 
element  method.  The  equations  have  been  solved  in  the  first  quadrant  shown  in  Figure  2(b) 
by  using  the  fundamental  solutions  given  by  eqns  (27a)  and  (27b).  It  has  been  assumed 
that  the  deviations  in  the  fields  of  stress  and  deformation  generated  by  the  interphase 
properties  of  the  dissimilar  fiber  0  extend  to  the  next  nearest  two  fiber-layers  only,  i.e.  to 
the  fibers  1  and  2  in  the  first  fiber-layer  and  to  the  fibers  3  -6  in  the  second  fiber-layer. 


Fig.  3.  Geometrical  construction  of  the  first  quadrant  fundamental  solutions;  ^  —  (si, ^2)  and 
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The  contours,  To,  r2,  Fj,  r4,  Fsand  F^,  contain  unknown  deviations  of  the  tractions 
and  displacements.  Along  these  contours,  the  integration  is  performed  clockwise  for  the 
matrix  side  integration  and  counterclockwise  for  the  fiber  side.  Suppose  that  the  contours 
Fo,  F,,  F2,  F3,  F4,  F5  and  F6  are  divided  into  Ao,  ^6  elements, 

respectively,  where  the  fields  are  taken  as  constants  over  each  element.  We  then  have  8A^o> 
8A^i,  8A^2,  8A^3,  8A^4,  8A^5  and  8A^6  unknowns  for  the  integrations  over  the  contours  Fq,  Fi, 
F2,  F3,  F4,  F5  and  Fa,  respectively  (two  traction  and  two  displacement  components  on  each 
side  of  a  contour).  Consequently  we  have  a  total  of  8(AoT  -f  A^3  + A^4-f  A^5  + A^6) 
unknowns.  The  same  number  of  equations  is  needed.  Equations  (21)  and  (24)  give 
4(A^o+A^i  + A^2  +  ^3  +  ^4  4"^5  +  ^6)  equations  and  the  continuity  of  tractions,  eqn  (9),  tog¬ 
ether  with  eqn  (13)  on  the  interface  yields  Ni~\- N2~\~ Ni-\- N4  +  Nsi-Nf)  equations. 

Hence  the  total  number  of  equations  is  the  same  as  the  total  number  of  unknowns.  Thus 
the  discretized  system  of  eqns  (21)  and  (24)  can  be  solved  numerically.  When  there  are  no 
interphase  flaws,  all  interphases  are  equally  divided,  specifically  in  1 20  or  60  elements,  such 
that  4No  =  =  2M2  ==  120  and  IN^  =  iV4  =  =  2^^  =  60,  When  Fo  contains  a  crack, 

the  crack  tip  element  and  the  element  ahead  of  the  crack  tip  element  are  further  divided 
into  smaller  elements  to  give  more  accurate  results.  The  number  of  elements  near  the  crack 
tip  is  increased  until  a  further  increase  does  not  change  the  numerical  results.  Numerical 
calculations  have  been  carried  out  for  FP/AI  composites  (see  Takahashi  and  Chou,  1988) 
with  the  following  properties : 

Aluminum  matrix  :  =  25.61  GPa,  =  0.345 


FP(Al203)  fibers:  /  =  1 57.9 GPa,  v‘  -  0.2. 

The  perfect  and  the  dissimilar  interphase  constants  {k,..  A,,)  iind  (A,.,  ko)  were  rendered 
dimensionless  by  division  by  where  a  is  the  radius  of  the  fibers.  In  the  computations, 
the  two  interphase  constants  were  taken  equal  in  magnitude,  and  thus  Ar/(/f“A0  ~  A^/ 
(/i"7^)  =  A.  Similarly  in  the  dissimilar  interphase  A,/(/f7^/)  ==  Ay(/f7^0  =  A’. 

In  the  initial  state  of  the  calculation  the  radial  interface  stress  a"’  (t[  —  t[  —  —/I") 

is  computed  under  the  assumption  that  eqn  (la)  applies  along  the  interface  of  the  dissimilar 
fiber.  If  a  negative  radial  stress  value  is  obtained  over  one  or  more  elements  and  the 
displacements  display  an  overlap,  eqn  (la)  is  replaced  by  eqn  (lb)  and  the  calculation  is 
redone  until  the  overlap  disappears  and  all  boundary  conditions  are  satisfied.  Values  of 
A  =  0.001,0.1,  I,  10,  lOOand  00  were  considered  for  the  interphase  stiffness  of  the  dissimilar 
fiber  and  k  =  00  for  the  neighboring  fibers.  The  fiber  volume  ratio  F,  was  chosen  as  0.2, 
0.4,  0.6  or  0.8. 


RESIJITS 

The  approximation  of  including  only  the  nearest  and  next-nearest  neighbors  in  the 
BEM  calculations  would  seem  to  he  a  reasonable  one.  The  question  does,  however,  arise 
whether  sufficiently  accurate  results  can  be  obtained  by  including  only  the  nearest  neighbors, 
i.e.  fibers  I  and  2  in  Fig.  2.  Intuitively  it  is  expected  that  this  simpler  approximation  will  be 
valid  for  smaller  volume  densities,  F,,  Ihe  fibers.  Test  calculations  have  been  carried  out 
for  the  case  of  uniaxial  tension,  for  Fj  =  0.2  and  \\  =  0,8,  and  A  =  10  and  A  =  I.  The 
results  are  shown  in  Fig.  4.  It  is  noted  that  the  inclusion  of  the  next-nearest  neighboring 
fibers  (solid  lines)  makes  a  small  diflerence  which  becomes  somewhat  more  pronounced  for 
smaller  values  of  A  and  larger  values  of  Fj.  The  numerical  results  presented  in  the  sequel 
have  been  carried  out  taking  into  account  only  the  next  nearest  neighbors. 

Figure  5  shows  the  stress  fields  tr^tro,  at  the  matrix  side  of  the  interphase 

and  the  normalized  strain  energy  density  U  in  the  interphase  of  the  dissimilar  fiber  for  the 
case  of  tension  in  the  .v  direction.  The  radial  stress  in  the  region  of  tensile  stresses  decreases 
as  the  interphase  stiffness,  A,  of  the  dissimilar  fiber  decreases.  The  maximum  value  of  the 
hoop  stress,  near  0  =  75'',  does,  however,  increase  as  A  decreases.  Other  calculations  not 
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Fig.  6.  ffy/ffo  at  the  matrix  side  of  the  interphases  of  the  dissimilar  fiber  0  and  the  neighboring  fibers 
1  and  2,  for  various  k  with  k  ~  oo  and  V{  =  0.6,  and  for  uniaxial  tension  in  the  x  direction. 


reported  here  show  that  the  magnitudes  of  (^rol^oy  and  0  at  the  dissimilar 
interphase  are  slightly  less  than  those  of  (^rol^oy  and  0  at  the  interphases  of  the 
perfect  composite  with  the  same  interphase  stiffness  (k  =  k)  in  all  interphases.  That  differ¬ 
ence  becomes  larger  for  lower  interphase  stiffness  k. 

Figure  6  shows  the  hoop  stresses  along  the  matrix-side  of  the  interphases  of  the  three 
fibers  0,  1  and  2  with  Ff  =  0.6  for  various^:  ^  =  oo,  10,  1, 0.1, 0.001.  Comparison  with  the 
results  of  Fig.  5  shows  that  the  hoop  stress  a^ja^  of  the  dissimilar  fiber  for  =  0.6  is  slightly 
smaller  than  that  for  Ff  =  0.4.  The  hoop  stresses  of  the  neighboring  fibers  1  and  2  slightly 
increase  as  k  decreases.  The  result  shows  that  for  low  k  the  positive  maximum  hoop  stress 
(^e)inax  occurs  near  6  =  70°  for  the  dissimilar  fiber,  near  0  =  —208°  for  the  neighboring 
fiber  1  and  near  0  =  —  30°  for  fiber  2. 

Figure  7  shows  the  hoop  stresses  along  the  matrix  side  of  the  interphase  of  the  dissimilar 
fiber  and  the  strain  energy  in  the  interphase  for  the  case  of  uniaxial  compression  in  the  x 
direction.  When  the  composite  is  compressed  in  the  x  direction  the  maximum  tensile  hoop 
stress  occurs  at  0  =  0°  for  low  k  and  at  0  =  90°  for  high  k  (above  10).  For  the  dissimilar 
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Fig.  7.  aJoQ  at  the  matrix  side  of  the  interphase  and  €  in  the  interphase  of  the  dissimilar  fiber  0,  for 
various  k  with  k  =  oo  and  V(  =  0.4,  for  uniaxial  compression  in  the  x  direction. 


fiber  the  maximum  strain  energy  occurs  near  0  =  55°  for  ^  =  1.  The  hoop  stresses  of 
the  neighboring  fibers  1  and  2  are  compressive  for  any  k. 

When  the  interphase  of  the  dissimilar  fiber  contains  interphase  cracks  as  shown  in  Fig. 
2(c),  the  behavior  of  ajao,  ObI^o  and  V  depends  greatly  upon  the  interphase  stiffness 
Ic.  Figure  8  shows  ag/ao  and  0  for  various  Ic  along  the  matrix  side  of  the  interphase  of  the 
dissimilar  fiber  for  an  interphase  crack  of  length  c  =  c/0.05236fl  =  9,  for  the  case  of  an 
uniaxial  tensile  stress  (To  in  the  x  direction.  The  maximum  hoop  stress  and  energy  density 
are  obtained  at  the  crack  tip.  The  magnitude  of  the  hoop  stress  increases  with  increasing 
values  of  k.  The  results  for  ag/ao  at  the  neighboring  fibers  1  and  2  are  also  shown  in  Fig.  8. 
The  overall  behavior  for  the  neighboring  fibers  is  very  similar  to  that  displayed  in  Fig.  6. 

Figure  9  shows  ffe/ffo  along  the  matrix  side  of  the  interphase  of  the  dissimilar  fiber  and 
0  in  the  interphase  for  ^  =  10  and  ^  =  1 ,  for  various  interphase  crack  lengths  e.  The  results 
show  that  for  *  =  10  the  positive  maximum  hoop  stress  occurs  at  the  crack  tip  and  increases 
as  c  increases.  On  the  other  hand  for  k  =  1  the  positive  maximum  hoop  stress  occurs  near 
0  =  75“  for  5  <  18.  Results  not  displayed  here  show  that  for  both  k  =  1  and  k  =  10  the 
stresses  ag  along  the  interphases  of  the  neighboring  fibers  1  and  2  increase  locally  as  C 
increases,  but  there  is  little  difference  between  k  =  1  and  k  =  10. 

A  stability  condition  is  needed  when  interphase  cracks  are  present.  The  strain  energy 
criterion  was  introduced  in  the  papers  by  Achenbach  and  Zhu  (1989,  1990)  to  investigate 
the  proclivity  of  an  interphase  crack  to  propagate.  Figure  9  shows  that  0  plotted  versus  6 
shows  maximums  at  the  tips  of  the  interphase  cracks  of  the  dissimilar  fiber.  For  ^  =  1,  the 
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Fig.  8.  at  the  matrix  sides  of  the  interphases  of  the  dissimilar  fiber  0  the  neighborinj  fibers 

1  and  2,  and  0  in  the  interphase  of  fiber  0,  for  various  k  with  k  —  oo  and  Vf  .  , 
uniaxial  tension  in  the  x  direction. 
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Fig.  9.  ,r.Kat  the  matrix  side  of  the  interphase  and  0  in  the  interphase  of  the  dissimilar  fiber  0  for 
various  £  with  fc  =  oo.  F,  =  0.4  and  £  =  10  and  £  =  1,  for  uniaxial  tension  in  the  x  direction. 
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maximum  values  first  increase  and  then  decrease  as  c  increases.  They  become  in  fact  smaller 
than  for  c  =  0,  and  hence  the  interface  flaw  will  be  arrested  for  that  case. 


CONCLUDING  COMMENTS 

A  fiber  whose  interphase  has  a  lower  stiffness  than  the  surrounding  fibers,  or  whose 
interphase  contains  cracks,  gives  rise  to  higher  fields  of  stress  and  deformation  in  its 
immediate  vicinity.  For  the  case  of  transverse  loading,  a  numerical  method  based  on  the 
solution  of  boundary  integral  equations  by  the  boundary  element  method  has  been 
developed  to  obtain  the  stresses  at  the  matrix  side  of  the  interphases  and  the  deformation 
energy  density  in  the  interphases,  for  the  dissimilar  fiber  and  its  nearest  neighbors.  The 
example  of  a  hexagonal  array  fiber  composite  has  been  discussed  in  some  detail,  and 
particular  attention  has  been  devoted  to  the  hoop  stresses  and  the  deformation  energy 
densities  since  critical  values  of  these  quantities  may  be  associated  with  radial  matrix 
cracking  and  interphase  cracking,  respectively. 
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APPENDIX  :  FIRST  QUADRANT  FUNDAMENTAL  SOLUTIONS 
The  components  of  the  two-dimensional  first  quadrant  fundamental  displacement  solutions  may  be  written 
as 


t/?.  (*, « )  =  f  A I  (*. «)  -  +  t/i  1  (*.  ?) 

=  +  -fej  J 


(At) 


Stress  states  at  neighboring  fibers 
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where  K  =  l/[87r^(l  —  v)].  Also 


=  Uu(x,i)-Vn(x,  =  ^1^; 


dR  dR  dRj  dR 
dxi  dx2  dxx  dx: 


!2  ^  ^  dR^~ 

2  dxx  dx2  dxi  dx2^ 


[dR  dR  dRj  dR2  dR^  dR^  dR,  dR, 

r/?, (X. ?)  =  t/2. (X. i)  +  u„ (X.  - ?) - 1/2, (X. -i)-u,dx,o  =  K  + 


y?2  (X,  {)  =  t/22 (X.  {)  +  f/jj (X,  -  I)  -  t/jj  (X,  -  «)  -  ^22 (X,  I) 


=  a:  1^(3- 


.  /f,«,  /5«V 

4v)ln— -~  +  L  )  +L^^.  )  L^J 


.  (A4) 


Here,  as  shown  in  Fig.  2,  R2,  and  R^  are  the  distances  between  the  load  point  in  the  first  quadrant  and  the 
image  points  in  the  second,  third  and  fourth  quadrants,  respectively.  Also,  R  is  the  distance  between  load  point. 
{,  and  the  field  point,  x.  Thus, 

/{  =  |x-{|,  f?2  =  l*+?l.  Ri  =  \x-l\.  (A5) 

The  first  quadrant  traction  solutions  may  be  expressed  as 

r?,  (X. «)  =  r, ,  (X,  {)  -  r, ,  (x,  - 1)  -  r,,  (x,  - « ) + r, ,  (x,  ?) 

where  C  =  —2ftK.  Also 

T?2  (X,  =  (X. «)  -  T,  2  (X,  - 1)  -  r,  2  (X,  -  {)  +  r,  2  (X,  ?) 


Cf^dRdRdR  „  .J  dR 

cr^dR2dR2dR2  .7  dRi 
R2  I  dn  dx,  5x2  ^ 

C  r  5«)  dffj  3«3  ,  dR,  a/JjNI 

~  R,\^  dn  dx,  dx,  ^’’^("'5x2  ”^5x,jJ 

C  [JR,  dR,  dR,  „  ^  J  dR,  dRA] 

+  T,  f  17r  ax,  ax2  ^  dx,  ax,  j J' 

r?,  (x.  {)  =  T2 ,  (X.  {)  +  r2 ,  (x,  - 1)  -  7-2 ,  (x.  -  ?)  -  Tj ,  (x,  ?) 


C  r  dR  dR  dR 
^  R  dn  dxi  dx2 


J  dR  dRV\ 


dn  dXi  dx2  V  ^^2 


-(l-2v)(«, 


C[^dRydR,dR,  ^  /  dR,  ^ 

'dx, 

R,  L  dn  dxi  dx2  \  ^^'2 
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and 


Here  the  subscript  /  denotes  the  load  direction  and  lakes  a  value  of  1  or  2,  while  the  second  subscript  indicates  the 
displacement  or  traction  direction. 


Cylinder  Theory  for  Unidirectional  Composite 

Materials 


Jyi-Jiin  Luo  and  Isaac  M.  Daniel 


The  theory  discussed  here  is  intended  to  provide  an  approximate  solution  of  the  mul¬ 
tiple  fracture  problems  of  unidirectional  fiber  reinforced  composite  materials.  The  strategy 
is  to  assume  a  reasonable  displacement  field  which  is  expected  to  be  good  enough  for  slowly 
varying  traction  boundary  conditions  on  the  lateral  surface  of  a  cylinder,  then  apply  the  vari¬ 
ational  principle  to  get  the  governing  equations  and  necessary  boundary  conditions.  All 
the  field  quantities(stresses,  strains,  and  displacements)  are  expressed  in  terms  of  average 
quantities  over  the  cross  section  of  the  cylinder.  The  unidirectional  multiple  fracture  prob¬ 
lem  of  a  unidirectional  composite  is  solved  by  patching  the  fiber  and  matrix  cylinders  to¬ 
gether. 


1  Basic  Theory 


For  a  hollow  cylinder  with  inner  radius  Ri  and  outer  radius  Ro,  the  displacement  field  is 
assumed  to  take  the  following  form 

Ur{r,z)  =  rfi{z)  + -f2{z), 

Uz{r,z)  =  w{z)  + g{r)K{z).  (2) 

The  radial  displacement  Ur{r,  z)  is  required  to  be  constant  at  r  =  Ro,  or 

Xlfi^Roy  —  ^o* 


If  we  set  the  displacement  at  r  =  Ri  to  u{z),  the  radial  displacement  Ur{r,  z)  can  be 

1 


written 

(3) 


1 


Note  that  the  radial  displacement  iirir,  z)  is  consistent  with  the  Lamd  solution  of  a  pressur¬ 
ized  cylinder. 

For  the  axial  displacement,  we  require  w(z)  to  be  equal  to  the  average  displacement  of 
the  cross  section  of  the  cylinder  A,  or 

g(r)dA  =  0,  (4) 

The  shear  strain  ')rz  at  the  outer  diameter  Ro  is  required  to  be  zero, 


'yrzi^J^Oi  2)  —  (^r,z  “I”  ^z,r)r^Ro  —  0* 

Substituting  the  displacement  fields  (eqs.  2  and  3)  into  the  above  equation,  we  have 


Equations  (4)  and  (5)  are  automatically  satisfied  if  g(r)  takes  the  following  form 


gif')  =  -  Rl\nr- a,  (6) 

where  a  is  defined  by 

a  =  jJfj-Rl\nr)dA,  (7) 

in  which  the  cross  section  area  of  the  cylinder  is  given  by  A  =  n\Rl  -  R]\. 

Putting  Ro  =  Rm  and  Ri  =  Rf  into  the  above  displacement  fields  for  matrix  cylinder 
of  a  unidirectional  composite  element,  we  get 


Urir,z) 
Uzir,  z) 


(-r  +  ^'^RjUmiz)  -f  -  ^^RmUo 


Rl 


R) 


•Wmiz)  -b  9mir)Kraiz), 


where  gmif')  is  given  by 


gmir)  =  -  Rl,\nr-  am, 

=  -  Rm^ffr)dA. 


(8) 

(9) 

(10) 

(11) 
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The  function  is  interpreted  as  the  radial  displacement  on  the  inner  surface  of  the  mat¬ 
rix.  Uo  is  the  radial  displacement  on  the  outer  surface.  The  function  Wm{r)  is  the  average 
axial  displacement,  and  Km{r)  is  interpreted  as  the  curvature  of  axial  displacement  at  large 
r. 

In  the  fiber  region,  set  Ro  =  0,  Ri  =  R/,  and  Uo  =  0,  then  eqs.  (2)  and  (3)  become 


Ur{r,  z) 


,(r,z)  =  wj{z)  +  gf(r)K{z), 


where  gf{r)  is  given  by 


9j{r)  = 

The  function  uf{z)  is  interpreted  as  the  radial  displacement  on  the  lateral  surface  of  the 
fiber.  The  function  Wf{r)  is  the  average  axial  displacement,  and  Kf{r)  is  interpreted  as  the 
curvature  of  the  axial  displacement  at  the  center  (r  =  0)  of  the  cross  section  of  the  fiber. 


1.1  Cylinder  Strains  and  Cylinder  Stress  Resultants 


The  strains  in  the  hollow  (matrix)  cylinder  can  be  found  from  the  displacement  fields  given 
by  eqs.  (2)  and  (3) 


Cr  =  Ur 


Cfl  = 


Rl-RI 


dw(z)  ,  .  dK(z) 

dg(i')  r  /  \  duiz) 

7r.  =  «r,.  +  U,,,  =  -^[K(z)-  ^Z 

Define  cylinder  strains  ci,  ej,  7,  </i,  and  92  by 


=  ^  I  t^dA, 

A  Ja 
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1  /  tr  +  eg  J  ^ 

-  aJa  2 

(21) 

(22) 

•  -  i /.¥'’•■"' 

(23) 

92  =  J^dir^zdA, 

(24) 

where  h  and  h  are  given  by 

/.=/,bwivA 

Substituting  eqs.  (16-19)  into  (20-24),  we  get  the  cylinder  strain-displacement  relations 

dw 

dz 

(25) 

^2  *—  j^2 

(26) 

4 

^  ^ /?2  RoU{z)\, 

-  Rf) 

(27) 

Ri  du{z) 

9i  -  «(^)  7^2  ’ 

(28) 

dK(z) 

®  -  ’ 

where  m  is  defined  by 

1  ^  RoRi  f  1 
m  A  Jav^ 

(29) 

Define  cylinder  stress  resultants  iVi ,  A^2)  -S',  Mi ,  and  M2  by 

ATi  =  ^  cr^dyl, 

(30) 

f  (Tr  +  (^e  j. 

N2  =  1  2 

(31) 

» -  /, 

(32) 

M.  = 

Ja  dr 

(33) 

M2  =  ^  g{i')(yzdA. 

(34) 
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Using  linear  elastic  stress-strain  relations  between  stresses  in  eqs.  (30-34)  and  strains  in  eqs. 
(20-24),  we  obtain 


N\  —  (‘^4"  2/x)/lci  4"  2Ai4c2, 

(35) 

N2  —  2(A  4"  4"  A^ci, 

(36) 

S  =  ^A-y, 

(37) 

Ml  =  ixhqi, 

(38) 

M2  =  (A  4-  2fj.)l2q2, 

(39) 

where  A  and  n  are  Lamp’s  constants. 

Substituting  the  cylinder  strain-displacement  relations  (eqs.  25-29)  into  the  cylinder  stress 
strain  relations  (eqs.  35-39),  we  get  the  cylinder  stress-displacement  relations 

^1  —  {X  +  2fl)A  ^ 

(40) 

(41) 

=  ^(^RiUo-Rou(z)y 

(42) 

M.  =  2*)’ 

(43) 

M,  =  (A  +  2^)//f>. 

(44) 

The  strain  energy  per  unit  length  of  the  cylinder  Uc  is  defined  as 

Uc  =  1  UdA 

JA 

(45) 

=  2  4-  Ctf  -1-  4-  ^)  dA. 

(46) 

Using  the  cylinder  strain  definitions  (eqs.  20-24)  and  the  cylinder  stress-strain  relations 
(eqs.  35-39),  (4  can  be  simplified  to 


Uc  =  —  |Ai4(ei  H- 2^2)^ -1- 2/Xv4(€j -j- 2^2) -f  -j- ^/i5j 

+  2fi)l2ql^  i  (47) 

=  —  -|-  2N2€2  4"  fn^S'y  -f  Miqi  M2q^ .  (48) 
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The  equations  discussed  in  this  section  can  be  applied  to  the  matrix  region  if  Ro  is  replaced 
by  Rm,  and  by  R /.  For  the  fiber  region,  Ro  and  Uo  are  replaced  by  zero,  and  Ri  replaced 
by  Rf. 


2  Principle  of  Virtual  Work 


From  eq.  (47),  we  have 

SUc  =  Ayl(ei  -|-  2c2){S(.i  -+■  25c2)  4"  -t-  462^^2) 

+m^fiA'yS'y  +  lihqiSqi  +  (A  +  2^)/2<]'2<^92-  (49) 

Using  the  cylinder  stress-strain  relations  (eqs.  35-39),  eq.  (49)  becomes 

SUc  =  NiSti  +  2N2Se2  +  m'^SS'^  +  MiSqi  4-  3/2^92-  (50) 


Now  consider  the  independent  virtual  displacements  Su{z),  5w{z),  Sk{z),  and  Suo  from  an 
equilibrium  state.  Using  the  cylinder  strain-displacement  relations  (eqs.  25-29),  the  vari¬ 
ation  of  strain  energy  SUc  can  be  expressed  as 


SUc  =  -^^{-‘^^2  +  ^-2mSRo/R,)Suiz)-^Sw{z) 

+  (m,  -  (iNMUi  +  2ms)iu, 

+  (51) 


The  virtual  work  SW  is 

SW  =  f  (TrSUr  +  TJu,)dS 
Js 


=  f  Tr{Ro,  z)Suo2it Rodz  +  f  Tr{Ri,z)Su{z)2nRidz 
Jzi  Jz\ 

4-  J  T^{Ri,z)[Sw{z)  +  g{Ri)SK{z)]2nRidz 

f  .  RiRo  \f  Rl\Su{zi)  (  R]\Suo]j. 

+ L  l(-^ + f)  4r  1 

,  fmf  \  \f  ,  RI\Su{z2)  ,  R]\SUo 

+ L  aftr?  Ir’^  V  J  -rt + -  V  J  -rt 


(52) 
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(53) 


+  /  T^{r,zi)[5w{zi)  +  g{r)5K{zi)]dA 
J  A 

+  /  T^{r,Z2)[Sw{z2)  +  gir)5K{z2)]dA. 

J  A 

The  principle  of  virtual  work  requires  that 

SU-5W=  r  SUcdz  -  8W  =  0.  (54) 


Since  Su{z),  5w{z),  8k{z)  are  arbitrary  in  the  domain  (01,^2).  fromeqs.  (51),  (53)  and  (54), 


we  conclude  that  the  equations  of  equilibrium  in  the  z^  <  z  <  Z2  domain  are 

-;r^;^(-2N2  +  ^-2mSR,/R.)  =  2l,  RiTr{Ri,  z),  (55) 

Rj-  Kf\  dz  J 

=  2nRiT,iRi,z),  (56) 

dz 

+  =  2nRig{Ri)URi,z).  (57) 

dz 

At  2:  =  21,  we  have 

''Mi{zi)  =  -f  ^^^Trir,zi)dA  or  (5tt(2i)  =  0,  (58) 

Ja  dr 

Niizi)  =  -f  T;,{r,zi)dA  or  ^«)(zi)  =  0,  (59) 

JA 

4^2(21)  =  -/  g{r)T^{r,zi)dA  or  4k(zi)  =  0.  (60) 

Ja 


The  boundary  conditions  at  22  are  similar.  From  8uo  we  have  the  constant  displacement 
condition 


J  ^N2  +  mSRi/Ro^dz  =  A  J  Tr(Ro,z)dz 

(r-^)(r,(r,2i)  +  r,(r,22))dA  or  =  0,  (61) 

which  implies  that  we  can  either  specify  iio  or  determine  Uo  from  given  tractions  on  the 
surface. 


Substituting  the  cylinder  stress-displacement  relations  (40-44)  into  eqs.  (55-57),  we  can 
get  the  equilibrium  equations  in  terms  of  displacements.  In  the  matrix  region,  these  equa¬ 
tions  are  , 


T^Rf^mhin  «  /  \  ,  \  /  /  \  H’mhm  t  ^  ^ 

- - “m(^)  +  >^mWjz)  -  K,„(2)  -  A,„  +  +  p2  ) 


2'nRj 


2A. 


Rl 


Ry\ 


tt, 
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(62) 


_  qn  /  r>  \  “1“  2^^^) 

—  -”irm(^/) -2^)  -  , 

/\m 

(Ayn  +  2/i  m)-4m  // /  >^  ^  //\  /n/T~>  \ 

2^  Jlj  ■”  Am^ynV'^)  ^  -2:), 

where  =  7r(/?^  —  /2y)  is  the  matrix  cylinder  cross-sectional  area. 
Similarly,  the  constant  displacement  condition  (eq.  61)  becomes 

2(A,n  -|-  ..  ,  \  ^771(^2)  —  w;m('2^l) 


22  -  2i 


R,(l  -  f){z2  -  z,)  i  + 

Am  +  /*m(l  +  /)  2no  ^  1  m  /  r>  -  \j. 

i_/-  p  /  T-rm\Rmi  Z)dz 

J-  J  -2^2  2^1 


i-m  -2^2  2^1 

RV 


Z2-Zi)  Ia^  ("'  /) 


2^m(^2  -  Zi) 

Applying  the  same  procedure  in  the  fiber  region,  we  have 


-  X,w',(z)  +  f^K',{z)  -  lihq>‘lK,(z)  =  -Z,(R„z), 


(Xf  +  2/jy)-^uiJ(z)  +  A/lly(z)  =  -T,j(Rf,z), 

R3 

(A/  -I-  2flj)-^K'l{z)  -  flfu'j(z)  -  flfR/Kfiz)  =  -TrfiRf,  z). 


The  above  equations  can  be  rewritten  as 


|<(fl  +  2(a„  -  l)<(e)  -  -  (2  +  ^)u„(() 


/ 

C2am(l  -  /)  /, 


(63) 

(64) 


(65) 

(66) 

(67) 

(68) 


R/Rrm(0 

'  Uo  \ 

Mm 

i-f  ( 

~Rm) 

(69) 

'!({)  -  2(a„ 

-iK(0 

_ 

RjTzmii) 

5 

(70) 

Mm 

,«)  +  <(() 

1 

1 

C2 

=  • 

RjTzmii) 

Mth. 

(71) 
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where 


f=RVRi, 


Al  ’ 


ASir^hm 
Al  ’ 


t  —  zj Rf^  ^m(0  o  A  ^rn{^)  ““  /^m (^/)^m (^)5 

Z/lm 

TrmiO  =  Trm{Rf,z),  T^O  =  nm{Rj,z). 
The  constant  displacement  condition  (eq.  65)  becomes 

4(0771  —  l)/_  ,  e^f  .  \  0^771(^2)  ~ 

— (wr“”  +  2(a„  -  1) - - = 


\1  f  /  ^Rm'  /^77i(s2  Cl)  *^€1 
r  Rf\ 


2/im(<f2  \Rf  f' 


Trm{r,  Zi)  +  Trmir,  Z2)  )dA, 


where 


“6-6  1. 


The  traction  terms  in  the  above  equation  are  zero  if  we  consider  the  symmetry  of  the  multiple 
crack  problem  of  unidirectional  composite  materials. 


Similarly,  in  the  fiber  region,  we  have 

ju'KO  -  2{a,  -  lyo'jiO  +  «.}({)  -  (4117  -  2)117(6  = 

07«,;(f)  +  2{a7-l)u'7(6  = 

-  “7(0  -  ‘l■>7(f )  = 


R,Tr,(() 

RfT.fiO 

RMO 


where 


TrfiO  =  Trj{Rj,  z),  T,fi0  =  nARf,  z). 
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3  Multiple  fracture  of  unidirectional  composite 

The  above  theory  can  be  applied  to  the  multiple  fracture  problem  of  unidirectional  compos¬ 
ites.  The  following  interfacial  conditions  are  discussed: 

1.  Perfect  bonding: 

Urf{Rf,  z)  =  Urm.{Rji  ^)) 

or 

Uf{z)  =  Um{z),  W/{z)  -h  Vf{z)  =  Wm{z)  +  Vm{z). 

2.  Constant  friction: 

Urf{Rf,  Z)  =  Urm{R/,  z),  T,j{Rf,z)  =  ±r„, 

or 

Uf{z)  =  Um{z),  Tzfii)  =  ±To. 

3.  Coulomb  friction: 

llr/{Rj,z)  =  Urm{Rj,  z),  Tzf{Rf,z)  =  ±nTrf{Rf,  z), 

or 

Uj{z)  =  Umiz),  ZfiO  =  ±fiTrf{0. 

4.  Noncontacting,  or  traction  free  interfaces: 

Trji^RmiZ)  =  0,  Trm{Rf^z)  —  0 


or 

TrfiO  =  0,  T,fiO  =  0. 


The  continuity  of  tractions  requires  that 


Trf{R/,z)  =  -Trm{Rf,z),  T2j{Rj,z)  —  -T^miRf^z), 


or 

TrfiO  —  TzfiO  —  ~TzmiO- 
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Consider  any  two  of  the  above  four  cases  connected  together  at  zi,  which  could  include 
fiber-matrix  debonding.  The  boundary  conditions  at  the  ends  zo  and  22  specify  either  the 
cylinder  displacements  or  the  cylinder  stress  resultants,  that  is,  at  ^  =  ^0  and  ^  =  6.  we 
have 


Uy 

or 

Mu 

specified, 

Um 

or 

Mim 

specified, 

Wf 

or 

Niy 

specified. 

Wm 

or 

Mm 

specified. 

Vj 

or 

A/2/ 

specified. 

Vm 

or 

M2m 

specified. 

The  continuity  conditions  at  zi  specify  the  continuity  of  cylinder  displacements  and  cylinder 
stress  resultants,  that  is, 

“/(6-)  =  W/(6  +  ),  ^/(^i-)  =  Wy(^i  +  ),  V/(^l-)  =  V/iCl-h), 

Miy(^l-)  =  Mim(^l+),  Niy(^l-)  =  Ntm(^l+),  M2/(6-)  =  ■^2m(6+)> 

Miy(^l-)  =  Aflm(6+))  ^1/(6“)  =  ^lmUl  +  ),  Af2/(6-)  =  Af2m(6+)- 


4  Crack  Opening  Displacements  and  Stress-Strain  Rela¬ 
tions 

Now  we  derive  an  important  result  from  the  cylinder  model.  Rewriting  the  equations  of 
equilibrium  (55-57)  in  the  fiber  region,  we  have 

=  2nR,T„(Ri,z),  (77) 

dz 

+  =  2nR,g,(R,)%,(R,,z),  (78) 
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and  in  the  matrix  region, 

Rl-R) 


dA/ im 
dz 


dM2m 


q  Rm  \ 

Rf) 

=  2nRfTrm 

(79) 

dN\m 

dz 

=  2TTRjTzm 

(80) 

■  +  Mim 

=  2tt  Rfgm{R])Tzm.{Rs,z). 

(81) 

The  constant  displacement  condition  (eq.  61)  gives 

J  ^N2m  I Rm^dz  =  Ajn  J  Trm{Rrrn^)dz 

+i  (r  -  S)  ^r,„(r,  z.)  +  T™(r, 


-  )  Zl)  +  Trmir,  22))  d  A. 


Now  we  focus  our  attention  to  the  global  behavior  of  the  representative  element  by  con¬ 
sidering  average  properties  over  [zi ,  Z2].  Integrating  eqs.  (77)  and  (80)  over  [zi,  Z2]  with  the 
weight  function  {z  —  zi),  one  gets 

r  Nif{z)dz  =  {z2-zi)Nifiz2)  +  2nRj  r  {z-zi)T,jiRf,z)dz  (83) 
Jzi 

f  NiTn{^)dz  =  {Z2  --  Zi)Nim{^2)  ~  ^l)Tzm{Rf -i  ^)dz,  (84) 

Jzi 

Integrating  eq.  (76)  over  [zi,Z2]y  one  gets 

/  ^  N2f{z)dz  =  l:Mif{z)  +  nR)  f  Trj{Rf,z)dz.  (85) 

Jzi  I  Zl 

Integrating  eq.  (79)  over  [21,22],  and  using  eq.  (82),  one  gets 

J  ^2miz)dz  =  2  Ja 

+  £’  lrR^iTr,^(R,,  z)  +  2)]  dz,  (86) 

r  S„(z)dz  =  [7’,„(r,20  +  r.„(r,2,)]irf>l 

_zR^  J»  ^  z)]&  (87) 
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Consider  now  a  representative  composite  element  with  matrix  cracks  at  z  —  zi  and  z  —  Z2. 
Due  to  symmetry,  the  following  boundary  conditions  are  required: 

Trf(r,Zi)  =  Trm(r,Z2)  =  0,  (88) 


and 

which  states  the  fact  that  there  is  no  global  stress  acting  on  the  lateral  surface  of  the  compos¬ 
ite.  Note  that  the  constant  displacement  at  the  matrix  outer  surface  is  automatically  satisfied 
in  that  displacement  assumption.  The  continuity  of  stresses  at  the  interface  requires 


I  ^  TrmiRm,z)dz  =  0, 

Jz\ 


(89) 


T,ARj,z)  =  -TURf,^\  (90) 

nj{Rf,z)  =  -T.m{Rf,z)-  (91) 


Debonding,  sliding  and/or  sticking  could  happen  at  the  interface.  These  conditions  usually 
would  require  further  elaboration,  but  intentionally  leaving  out  these  conditions  would  im¬ 
ply  that  the  results  are  applicable  to  all  of  them.  Applying  the  above  boundary  conditions 


(eqs.  90  and  91)  to  eqs.  (83-87),  one  gets 


1 

f  Nij{z)dz 
Jz\ 

Z2-  Zi  . 

1 

1  Nim{z)dz 

Z2  —  Zi  J 

1 

f  N2j{z)dz 

Jzi 

Z2-  Zi. 

1 

f  N2m{z)dz 

Z2  —  Zi  J 

1 

/  Sm{z)dz 
Jz\ 

Z2  -  Zi 

Ni/{z2)  +  AjQ, 

(92) 

Nlm{z2)  —  AjQi 

(93) 

AfP, 

(94) 

-AjP, 

(95) 

TT  Rm  Rf  p 
m 

(96) 

where  the  average  interface  radial  traction  P,  and  the  normalized  first  moment  of  interface 
longitudinal  traction  Q  are  given  by 


P  = 


Q  = 


1 


Z2  -  Zi 
2 

Rf  Jz\ 


r  Trj{Rj,z)dz  = - —  r  Trm{Rf,z)dz, 

Jz\  ^2  ^ 

Z2  —  Zi  tij  Jzi  Z2  —  Zi 


(97) 

(98) 
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Our  intention  is  to  derive  average  stress  and  strain  relations  by  integrating  field  quant¬ 
ities  over  the  ^-axis.  The  average  strains  are  calculated  by  integrating  the  cylinder  stress- 
strain  relations  over  [2:1,  Z2],  for  example, 


— —  /  eijdz 
Z2  —  Zi  Jz\ 

r  -  2u,N^j)dz 

Z2  —  Zi  Jzi  rjjAf 


Substituting  eqs.  (92  and  94)  into  the  above  equation,  one  gets 

p  +  g. 


Ejcif  = 


Following  the  same  procedure,  one  gets  the  other  average  strain  components: 

n  ^  ,  /I  _  P  _  ,,,n 


Efe2f  =  -^j 


+  (i  ~  ~  ^jQi 


Em^ltn  — 


Em^2r> 


^  Af  '  ^ 

Nirn{^2)  ,  ‘^f^m  jj  ^ 

w  i-r’ 

Nxm{z2)  fil- I'm)  a  ,  ^ 

-‘'’'-j: - 

TT  Rm  Rf  rj 


(100) 


(101) 


(102) 


(103) 


where  all  quantities  with  overbar  represent  their  averages  over  [zi,  Z2]. 

When  the  matrix  has  lower  strength  and  suitable  stiffness,  the  matrix  cracks  will  develop 
first,  and  parallel  matrix  cracks  could  occur.  Consider  the  problem  with  matrix  cracks  at  zi 
and  Z2,  see  fig.  1.  The  interface  may  have  debonding  and/or  sliding.  The  composite  longit¬ 
udinal  strain  ei  is  the  same  as  the  average  fiber  longitudinal  strain  hj,  or,  using  eq.  (99), 


^1/(^2)  Ip 


(104) 


The  composite  transverse  strain  €2  is  given  by 


(105) 
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Integrating  the  cylinder  strains-displacement  relations  (26)  and  (27)  for  the  matrix  cylinder, 
the  above  equation  can  be  rewritten  as 


^2  —  ^2m  — 


mfAr, 


-In 


Substituting  eqs.  (102)  and  (103)  into  the  above  equation,  we  have 


_  _  I'mNimjz-i) 


The  average  matrix  crack  opening  displacement  (X)D„  is  defined  by 

00D„,  =  [Wj{z2)  -  Wj{zi)]  -  [tU,n{z2)  “ 

=  (^2-  Zi){Cif  -lim). 


(106) 


(107) 


We  define  the  number  of  matrix  cracks  per  unit  length,  or  the  matrix  crack  density,  A,„  by 


Am  — 


(108) 


Combining  eqs.  (99),  (101)  and  (108),  one  can  express  the  normalized  matrix  crack  density 
AmCDDm  by 


AmOODm  —  Cl/  —  fli 
^  Nlf{z2) 

EjAj 


+ 


N\m{z2) 

fp  A 

fl^m 


lEj  il-f)E„ 


2P  + 


+ 


f 


lEf  {l-f)E^\ 


Q.  (109) 


The  average  gap  between  the  fiber  surface  and  the  inner  matrix  surface  is  given  by 


—  2Amfc/C0D,', 


where  the  average  displacement  u/  and  are  given  by 

_  1 

uf  =  - -  /  uj{z)dz, 

Z2  —  Zi  Jzi 

_  1  fZ2 

Urn  =  -  /  Um{z)dz, 

Z2  —  Zi  Jzi 


(110) 


(111) 

(112) 
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OOD,  is  the  average  interface  crack  opening  in  the  debonded  region,  and  id  is  the  debonded 
length. 


Integrating  the  strain-displacement  relations  (25—27)  for  both  the  fiber  and  matrix  re¬ 
gion,  one  can  rewrite  the  average  gap  in  eq.  (110)  as 


2\midGODi  - 

Using  eqs.  (100),  (102)  and  (103),  the  above  equation  can  be  written  as 


(113) 


2A,„400D.-  _  fmNimjzT) 

■ 


EfAf 


RmAr 


1-i^f  i  +  f  +  {i-fh 
Ej  (!-/)£ 


'm^m 

'm 


P  + 


m 


Ej"  (\-f)E^\ 
Rearranging  eqs.  (104),  (107),  (109)  and  (114)  in  matrix  form,  we  have 


Q.  (114) 


Cl 

^2 


}-\i 


Ef  0 


0 

"  Em  -■ 


Nh(z2) 

Am 


1 


-t- 


.Hii 

Ef 

AL 


'(l-l)Em  (l-/)Em  J 


i;i' 


(115) 


and 


I 


A„CDD,„ 

2A,„fdOOD./i?/ 


1 

Ef 

}iL 


E 


_ l_ 

Em 

_ 

Em 


Nxf(z2) 

Ef 

I^lm{z2) 
Em 


+ 


_  ji'm  1  _| _ L _ 

E,  (\-J)Em  Ef  ^  (l-f)Em 

l—ff  _  l-i-/+(l~/)t'm  I  ft'm 

■  Ef  (\-J)Em  Ef  (l-})Em- 


(116) 


Eliminating  P  and  Q  from  eqs.  (115)  and  (116),  after  complex  algebraic  operations,  one 
gets  the  following  results 


{ 


(1  - 

2XmidOODi/Rf 


1  _  &11  bi2  f 
J  ^21  ^22  \ 


Cl  ~  ■p" 
C2  + 


where 


^  J  /(t^/  -  t'm)  ^  /(I  +  t^/)(l  -  I  /i/ 


bi2  — 


1  -  '  1  ~ 
(1  -  f)i^J  -  ^rn) 


1  - 


(117) 

(118) 
(119) 


m 
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(120) 


&21 

&22 


fit's  -  I'm)  _  ft'mjl+l'f)  / if/  _  A 
1  -  l/fl/m  2(1  -  l/fUm)  \fim  ) 
^  -  Djt'f  -  I'm) 

1  -  I'ft'm 

(1  -/)(!  +t^>n)(l  -2t//)  A  _ 

2(l-J//l/^)  V  fif  ) 


(121) 


The  parameters  bij  (i,  j  =  1 , 2)  are  the  elastic  properties  of  the  composites,  and  depend  only 
on  the  constituent  elastic  constants  and  fiber  volume  ratio;  ai  is  the  composite  longitudinal 
stress  given  by 


_  _  +  -^Im 

+  /Im 


(122) 


Note  that  iVim(^2)  is  zero  because  of  the  presence  of  matrix  cracks;  Ei  is  the  undamaged 
composite  Young’s  modulus  given  by 


El  —  f  Ef  -j-  (1  —  f)Em  + 


_ 4/(1  -  f)it'}  -  UmfUm _ 

1  +  /(I  -  2t/m)  +  (1  -  /)(!  -  2uj)n,n/liJ  ’ 


and  1^12  is  Poisson’s  ratio  given  by 


1^12  =  ft's  +  (1  -  /)^'m  + 


/(I  -  f)i^S  -  ^m)[l  -  2t/,n  -  (1  -  2vj)Hmly.s\ 
1  +  /(I  -  2vm)  +  (1  -  /)(!  -  2uj)nmlsts 


(123) 


(124) 


5  Results  and  Discussion 


Introducing  the  experimentally  determined  constituent  and  composite  properties  in  eq.  ( 1 1 7), 
we  obtain  results  for  the  total  crack  opening  per  unit  length  COD™  and  for  2Am  AOOD,/ Rj, 
which  are  plotted  versus  applied  strain  ei  in  Fig.  2  and  3.  which  are  the  results  from  two  dif¬ 
ferent  batches  of  SiC/CAS.  They  both  show  that  interfacial  debonding  happens  a  little  bit 
after  matrix  cracks  start.  Fiber  breaks  are  ignored  in  the  calculation.  Since  fibers  only  break 
under  higher  loading  (mostly  after  matrix  crack  saturation),  these  results  are  believed  to  be 
acceptible  for  lower  strains,  and  can  be  used  only  for  reference  under  higher  strain. 
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Fig.  1:  Unit  Cell  Between  Matrix  Cracks  I 
00D,„  and  Interface  Debonding  Cr 


Stress,  a-\,  (ksi) 


Strains,  ei,  €2,  (%) 


Fig.  2:  Stress  Strain  Curves  and  Predicted  Crack  Opening  Displacements  OOD^  and  COD, 
of  Unidirectional  SiC/CAS  Composite  under  Longitudinal  Tension. 
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CODm,  2A^<,  COD/R,.  (%) 


Strains,  e^,  62,  (%) 


Fig.  3:  Stress  Strain  Curves  and  Predicted  Crack  Opening  Displacements  000^  and  OOD, 
of  Unidirectional  SiC/CAS  Composite  under  Longitudinal  Tension. 
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COD^,  2An/a  COD/R,,  (%) 
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Abstract— Failure  mechanisms  were  studied  under  the  microscope  in  a  crossply  silicon  carbide/glass- 
ceramic  composite  under  axial  tensile  loading.  Failure  initiation  takes  place  in  the  90°  layer.  It  takes 
the  form  of  radial  matrix  cracks  around  the  fibers,  followed  by  interfacial  cracks,  which  in  turn 
coalesce  into  transverse  macrocracks.  These  transverse  macrocracks  in  the  90°  layer  reach  a  charac¬ 
teristic  saturation  crack  density  with  a  minimum  crack  spacing  of  the  order  of  the  layer  thickness. 
Subsequently,  transverse  matrix  cracks  are  generated  in  the  0°  layer,  increasing  in  density  up  to  a 
minimum  crack  spacing  of  the  order  of  eight  fiber  diameters.  This  stage  of  failure  is  accompanied 
by  fiber-matrix  debonding  and  some  fiber-failures  in  the  0°  layer.  In  the  third  stage  of  damage 
development,  the  macrocracks  of  the  90°  layer  branch  off  and  connect  with  the  0°  layer  cracks  in  a 
characteristic  “delta”  pattern.  This  is  finally  followed  by  delamination  and  additional  cracking  in 
the  90  layer  prior  to  ultimate  failure.  The  various  lailure  mechanisms  and  their  interactions 
were  discussed  and  compared  with  predictions  of  prior  experimental  and  analytical  studies  of 
unidirectional  and  crossply  composites. 


1.  INTRODUCTION 

The  newer  generation  of  advanced  composites  designed  to  meet  high  temperature  per¬ 
formance  requirements  includes  ceramic-matrix  composites.  Ceramic  composites  display 
appreciable  ductility  and  have  higher  fracture  toughness  and  better  fatigue  endurance  than 
comparable  monolithic  ceramics.  Promising  reinforcements  include  carbon,  silicon  carbide 
(SiC)  and  mullite  fibers.  Matrices  used  successfully  include  glass,  glass-ceramic,  carbon, 
silicon  carbide,  silicon  nitride  and  aluminum  oxide.  Composites  consisting  of  glass-ceramic 
mauices,  such  as  lithium  aluminosilicate  (LAS),  calcium  aluminosilicate  (CAS)  and  mag¬ 
nesium  aluminosilicate  (MAS)  reinforced  with  silicon  carbide  yarn  or  monofilament  have 
been  developed;  see  Prewo  and  Brennan  (1982),  Brennan  and  Prewo  (1982),  Larsen  and 
Adams  (1989).  These  materials  can  be  fabricated  in  both  unidirectional  and  multidirectional 
laminate  form,  and  can  attain  strengths  of  1000  MPa  (140  ksi)  and  ultimate  strains  of  over 
1%.  Fracture  toughness,  as  determined  by  a  notched  beam  method,  is  more  than  five  times 
that  of  the  monolithic  ceramic,  Brennan  and  Prewo  (1982). 

The  micromechanics  of  stress  transfer  and  fracture  have  been  studied  experimentally 
and  analytically  by  many  investigators  for  the  case  of  unidirectional  composites,  starting 
with  the  pioneering  work  by  Aveston  et  al.  (1971)  and  Aveston  and  Kelly  (1973).  Fun¬ 
damental  approaches  used  include  shear  lag  or  stress  transfer  models  and  fracture  mechanics 
with  stress  or  energy  fracture  criteria.  Failure  mechanisms  in  a  SiC/LAS  unidirectional 
composite  under  longitudinal  loading  have  been  observed  and  analyzed  by  Marshall  and 
Evans  (1985,  1988)  and  Marshall  et  al.  (1985).  Further  analytical  studies  were  described 
by  Budiansky  et  al.  (1986),  Ilsueh  (1988),  Charalambides  and  Evans  (1989),  McCartney 
(1989),  Hutchinson  and  Jenson  (1990),  Zok  and  Spearing  (1992)  and  Weitsman  and  Zhu 
(1993).  The  authors  have  reported  both  experimental  and  analytical  investigations  of  a 
longitudinally  and  transversely  loaded  unidirectional  SiC/CAS  composite ;  see  Daniel  et  al. 
(1989a,b,  1992,  1993)  and  Lee  and  Daniel  (1992).  Most  studies  to  date  on  the  behavior  of 
composites  under  transverse  loading  are  analytical,  such  as  those  by  Hashin  (1983),  Benven- 
iste  (1985)  and  Achenbach  and  Zhu  (1990). 
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The  problem  of  crossply  laminates  has  been  studied  extensively  for  polymer-matrix 
composites.  The  primary  failure  mechanism  in  these  composites  is  transverse  matrix  crack¬ 
ing  in  the  90°  layer  reaching  a  limiting  crack  density.  This  is  followed  by  longitudinal  matrix 
cracking  in  the  0°  layers,  local  delaminations  at  the  intersections  of  matrix  cracks  and, 
finally,  fiber  fractures  in  the  0°  layer  leading  to  ultimate  failure,  Daniel  et  aL  (1988).  In  the 
case  of  brittle-matrix  composites  the  type,  sequence  and  interaction  of  failure  mechanisms 
are  different  and  characteristic  of  the  constituent  material  properties.  Some  experimental 
work  on  crossply  ceramic-matrix  composites  has  been  reported  by  Sbaizero  and  Evans 
(1986).  They  tested  symmetric  [0/90],  SiC/LAS  laminates  in  tension  and  bending.  However, 
the  observed  failure  modes  were  not  easily  explained.  The  initial  delamination  failure 
observed  would  only  be  possible  near  free  edges  in  the  case  of  a  material  with  very  low 
interlaminar  shear  strength.  The  orientation  of  the  transverse  ply  cracks  could  only  be 
explained  by  assuming  compressive  residual  stresses  in  the  90°  layers.  Subsequent  work  on 
crossply  laminates  dealt  with  a  fracture  mechanics  model  of  a  delamination  crack  formed 
at  the  tip  of  a  notch  cutting  through  0°  and  90°  layers,  see  Sbaizero  et  al.  (1990)  and 
Charalambides  (1991).  Mall  and  Kim  (1992)  also  conducted  experimental  investigations 
of  failure  mechanisms  in  SiC/CAS  laminates.  They  observed  that  initial  failure  consisted 
of  matrix  microcracks  occurring  and  growing  in  all  plies  in  a  random  manner. 

This  paper  deals  with  a  systematic  investigation  of  the  failure  mechanisms,  their  growth 
and  interactions  in  a  crossply  SiC/CAS  laminate  under  uniaxial  tension.  It  represents  the 
first  systematic  experimental/analytical  study  of  a  crossply  laminate  with  brittle  matrix. 
Previous  experimental  work  mentioned  before  reported  random  or  other  edge-dominated 
failure  mechanisms  with  no  clear  trends.  The  failure  mechanisms  observed  and  recorded  in 
the  present  study  are  substantially  different  from  those  reported  before  and  show  a  clear 
and  systematic  trend.  The  initial  stages  of  failure  in  each  layer  are  interpreted  on  the 
basis  of  the  failure  micromechanics  of  the  unidirectional  material  under  longitudinal  and 
transverse  loading.  The  interaction  of  failure  mechanisms  in  the  0°  and  90°  layers  and  the 
later  stages  of  damage  are  explained  on  the  basis  of  classical  lamination  theory  as  well  as 
on  damage  models  available  for  crossply  laminates.  The  various  failure  mechanisms  and 
damage  states  were  correlated  with  the  overall  mechanical  behavior  of  the  laminate.  The 
main  contribution  of  this  paper  lies  in  the  integration  of  the  micromechanics  of  brittle 
matrix  single  layers  and  the  macromechanics  of  a  crossply  laminate  in  predicting  and 
interpreting  the  type  and  sequence  of  failure  mechanisms  and  their  interaction. 


2.  EXPERIMENTAL  PROCEDURE 

The  material  investigated  was  SiC/CAS,  calcium  aluminosilicate  glass-ceramic 
reinforced  with  silicon  carbide  fibers,  manufactured  by  Corning  Glass  Works.  The  fiber  is 
silicon  carbide  yarn  known  as  Nicalon  (Nippon  Carbon  Co.).  This  fiber  is  available  in 
continuous  length  tows  with  an  average  diameter  of  1 5  /im.  The  composite  material  for  this 
study  was  obtained  in  the  form  of  unidirectional  and  [0/9()]2s  [O/9O2],  crossply  laminates. 

The  unidirectional  material  was  characterized  first  to  obtain  average  physical  and 
mechanical  properties.  Unidirectional  specimens  instrumented  with  strain  gages  were  tested 
under  longitudinal  and  transverse  tensile  loading  in  a  servohydraulic  testing  machine. 

The  macroscopic  response  of  the  crossply  laminates  was  measured  by  testing  specimens 
under  axial  tension.  The  specimens  were  15.2  cm  (6  in.)  long  and  1.27  cm  (0.5  in.)  wide. 
They  were  tabbed  with  3.81  cm  (1.5  in.)  long  crossply  glass/epoxy  tabs.  These  specimens 
were  tested  in  a  servohydraulic  testing  machine  (Instron)  using  an  extensometer  to  monitor 
strain  to  failure. 

Failure  mechanisms  were  monitored  by  testing  smaller  crossply  coupons  under  the 
microscope  with  a  specially  designed  fixture.  The  specimens  were  76  mm  (3  in.)  long  and 
12.7  mm  (0.5  in.)  wide  coupons.  Potting  epoxy  was  cast  around  these  coupons  to  help  in 
subsequent  polishing  and  loading  of  the  specimen.  One  edge  of  the  specimen  was  polished 
for  microscopic  observations. 
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Fig.  1 .  Setup  for  microscopic  observation  of  failure  mechanisms  in  composite  specimens  under  load. 

A  special  loading  fixture  was  designed  and  built  for  loading  the  specimens  under  the 
microscope  (Fig.  1) ;  see  Daniel  et  al.  (1989a,b).  Load  was  applied  and  controlled  by  means 
of  a  pneumatic  cylinder.  The  specimen  with  the  grips  attached  was  mounted  onto  a  reaction 
frame  attached  to  one  end  of  the  pneumatic  cylinder.  One  end  of  the  specimen  was  connected 
to  the  moving  piston  while  the  other  end  was  reacted  at  the  other  end  of  the  reacUon  frame 
through  a  strain  gage  load  cell.  The  entire  assembly,  including  air  cylinder,  reaction  frame, 
specimen  with  grips  and  load  cell,  was  suspended  with  a  counterweight  from  a  movable 
upright  frame.  Thus,  any  part  of  the  loaded  specimen  could  be  moved  to  the  stage  of  the 

microscope  without  any  weight  or  force  exerted  on  the  microscope.  .  i 

Two  types  of  tests  were  conducted,  continuous  and  intermittent  loading  tests,  n  e 
first  type,  loading  is  applied  continuously  to  failure.  A  specific  area  of  the  specimen  is 
observed  through  tlie  microscope  and  recorded  by  a  video  camera.  This  type  of  loading 
yields  a  continuous  stress-strain  curve  not  affected  by  stress  relaxation,  as  is  the  case  during 
intermittent  loading.  However,  only  one  small  area  of  the  specimen  is  being  monitored  and 
many  interesting  phenomena  occurring  outside  the  field  of  view  are  missed.  .  . 

In  the  second  type  of  test,  loading  is  applied  in  steps.  The  stress  and  strain  m  the 
specimen  are  recorded  in  the  beginning  and  at  the  end  of  each  loading  step.  Photo¬ 
micrographs  are  taken  at  difl'erent  areas  of  the  specimen  while  the  strain  is  held  constant. 
One  problem  in  this  type  of  testing  is  stress  relaxation  during  the  holding  period,  which 
becomes  more  pronounced  as  damage  progresses. 


3.  RESULTS 

Material  characterization  .  .  .  c  r, 

Some  basic  properties  of  the  matrix  and  fiber  constituents  obtained  from  the  literature 

and  from  previous  tests  are  tabulated  in  Fable  1. 

Table  1.  Constituent  material  properties 
Property  CAS  matrix!  liberj 

Maximum  use 

temperature,  °C  CF) 

Fiber  diameter  (yan) 

Density  (g  cm“^) 

Coefficient  of  thermal 

expansion,  I0"^‘^C“'  (10  ‘) 

Elastic  modulus, 

GPa  (10^  psi) 

Tensile  strength, 

MPa  (ksi) 

t  Larsen  and  Adams  (1989). 

JPrewo  and  Brennan  (1982),  Daniel  et  al  (1989). 


1350  (2460) 

1300  (2370) 

_ 

15 

2.8 

2.6 

5.0  (2.8) 

3.2(1.8)-4.9  (2.7) 

98  (14.2) 

170  (24.6) 

124(18) 

(flexural) 

1930  (280) 
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Strain,  c,  (%) 


Fig.  2.  Stress-strain  curves  to  failure  of  unidirectional  SiC/CAS  composite  under  longitudinal  and 

transverse  tension. 


Examination  of  photomicrographs  of  transverse  cross-sections  showed  that  the  fiber 
diameter  ranges  between  8  and  20  /tm,  with  a  median  value  of  15  /tm.  The  fiber  distribution 
is  nonuniform  with  an  average  fiber  volume  ratio  of  0.39.  A  low  degree  of  porosity 
was  measured  (Kv=  0.01). 

Typical  stress-strain  curves  to  failure  of  the  unidirectional  composite  under  longi¬ 
tudinal  and  transverse  tension  are  shown  in  Fig.  2.  Basic  properties  obtained  from  such 
macromechanical  tests  of  the  unidirectional  SiC/CAS  composite  are  tabulated  in  Table  2. 

Damage  initiation  ami  development  in  90°  layer 

The  first  stage  of  damage  development,  consisting  of  microcracking  in  the  90°  layers, 
was  studied  by  testing  [0/90]2s  crossply  laminates  under  the  microscope ;  see  Daniel  et  al. 
(1989a, b). 'The  first  microcracks  originated  at  the  fiber-matrix  interface  and  were  nearly 
normal  to  the  interface,  which  indicates  that  failure  is  caused  by  the  circumferential  tensile 
stress  in  the  matrix.  The  randomness  of  the  fiber  packing  makes  it  difficult  to  draw 
conclusions.  Flowever,  some  general  patterns  can  be  identified. 

When  the  fibers  are  closely  packed,  usually  in  a  near  hexagonal  array,  radial  cracks 
initiate  at  approximately  45°  from  the  loading  axis.  An  analysis  for  a  closely  packed 
hexagonal  array  of  fibers  showed  that  for  a  low  stiffness  interphase,  the  peak  circumferential 
stress  is  the  critical  one  and  it  occurs  near  the  45°  location,  Achenbach  and  Zhu  (1990). 
This  of  course  leads  to  radial  cracking  at  this  location,  as  observed  experimentally. 

When  fibers  are  farther  apart  and  are  surrounded  by  a  relatively  large  volume  of 
matrix,  radial  cracks  occur  at  approximately  90°  from  the  loading  axis.  For  a  low  stiffness 


Table  2.  Measured  properties  of  SiC/CAS  unidirectional  composite 


Property 

Value 

Fiber  volume  ratio,  Vf 

0.39 

Ply  thickness,  t,  mm  (in.) 

0.38  (0.015) 

Longitudinal  modulus,  GPa  (Msi) 

127(18.4) 

Transverse  modulus,  £"2,  GPa  (Msi) 

112(16.2) 

In-plane  shear  modulus,  (7,2,  GPa  (Msi) 

52  (7.5) 

Out-of-plane  shear  modulus,  C23,  GPa  (Msi) 

32  (4.7) 

Major  Poisson’s  ratio,  Vj2 

0.18 

Longitudinal  tensile  strength,  Fi„  MPa  (ksi) 

448  (65) 

Transverse  tensile  strength,  F21.  (^si) 

28  (4) 

Longitudinal  ultimate  tensile  strain,  c“, 

9x  10-^ 

Transverse  ultimate  tensile  strain,  £2, 

2.5x10“" 

Longitudinal  coefiicient  of  thermal  expansion, 

(/4e"F  '‘) 

4.05  (2.25) 

Transverse  coefficient  of  thermal  expansion, 
a2,  (/ic^’F"') 

4.23  (2.35) 

Out-of-plane  coefficient  of  thermal  expansion, 
o£3,  /ifi'K:-'  (/^e“F”‘) 

4.32  (2.40) 
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(C) 


Fig.  3,  Development  of  failure  mechanisms  in  transversely  loaded  ceramic  matrix  composite, 
(a)  Initial  radial  cracks  around  closely-packed  fibers,  (b)  Initial  radial  cracks  around  isolated  fibers, 
(c)  Interfacial  cracks,  (d)  Interconnection  of  radial  and  interfacial  cracks  and  formation  of  macrocracks. 


interphase  the  circumferential  stress  at  the  90^"  location  becomes  critical,  resulting  in  the 
observed  radial  cracks,  Daniel  et  aL  (1989). 

As  the  load  increases,  a  new  failure  mechanism  develops,  consisting  of  interface  cracks 
in  the  area  along  the  loading  axis  over  an  arc  20  about  the  0°  point  with  6  <  45°.  These 
interface  cracks  are  not  immediately  connected  with  the  radial  cracks  developed  earlier. 
Additional  radial  cracks  may  develop  at  this  stage  of  damage  development,  and  then,  radial 
and  interfacial  cracks  are  connected  to  form  a  long  continuous  crack  in  the  transverse  layer. 
The  various  failure  mechanisms  and  their  sequence  are  illustrated  schematically  in  Fig.  3. 

The  influence  of  these  failure  mechanisms  on  the  stress-strain  behavior  of  the  [0/90]2s 
crossply  laminate  is  shown  in  Fig.  4.  The  laminate  behaves  linearly  up  to  an  applied  stress 
of  27.6  MPa  (4  ksi)  when  the  first  radial  cracks  appear  in  the  matrix.  This  causes  a  sudden 
drop  in  stiffness  (secant  modulus).  At  a  higher  applied  stress,  32.4  MPa  (4.7  ksi),  interfacial 
cracking  occurs  causing  another  sudden,  but  smaller  drop  in  stiffness.  The  interconnection 
of  radial  and  interfacial  cracks  starts  at  a  higher  stress  (36.6  MPa ;  5.3  ksi)  and  progresses 
more  gradually  to  the  formation  of  continuous  macrocracks,  resulting  in  similarly  gradual 
reduction  in  laminate  stiffness.  Similar  stiffness  reductions  have  been  predicted  by  Zhu  and 
Achenbach  (1991)  for  a  transversely  loaded  composite  of  the  same  properties,  undergoing 
the  same  failure  process. 

The  stiffness  reductions  illustrated  in  Fig.  4  are  more  pronounced  in  this  test  due  to 
creep  effects.  The  loading  was  applied  intermittently  over  a  relatively  long  period  (more 


Fig.  4.  Stress-strain  curve  of  [0/90]2s  crossply  laminate  with  corresponding  stages  of  damage  devel¬ 
opment  in  transverse  layer. 
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Fig.  5.  Stress-strain  curve  of  [0/90^]^  SiC/CAS  crossply  laminate  under  uniaxial  tension. 


than  30  min)  to  allow  for  photomicrographic  recording  of  the  various  stages  of  damage 
development.  The  first  continuous  (across  the  90°  layer)  macrocracks  were  formed  at  an 
applied  stress  of  approximately  40  MPa  (5.8  ksi),  which  corresponds  to  an  axial  stress  of 
c72  =  38  MPa  (5.5  ksi)  in  the  90°  layer.  The  in  situ  lamina  strength  is  comparable  to  the 
value  measured  by  Mall  and  Kim  (1992)  and  is  higher  than  the  transverse  tensile  strength 
/'2t=  28  (4  ksi)  measured  from  a  90°  transverse  tensile  test.  Continued  loading  increases  the 
number  of  macrocracks  up  to  a  saturation  density,  corresponding  to  a  minimum  macrocrack 
spacing  of  the  order  of  the  90°  layer  thickness. 

Figure  5  shows  a  stress-strain  curve  for  a  [0/902]^  crossply  laminate  under  uniaxial 
tension.  Microcracking  in  the  90°  layer  (layer  2)  started  at  an  applied  strain  of  approxi¬ 
mately  0.025%.  Transverse  macrocracks  reached  a  minimum  crack  spacing  of  0.65  mm 
(0.0255  in.),  comparable  to  the  layer  thickness  of  0.68  mm  (0.027  in.),  at  an  applied  stress 
of  approximately  69  MPa  (10  ksi),  corresponding  to  a  strain  of  approximately  0.12%. 

The  initial  modulus  in  the  linear  region  is 

E,  =  117GPa(16.9Msi), 

which  agrees  with  the  calculated  laminate  modulus.  The  reduced  secant  modulus  at  the 
point  of  crack  saturation  (a^=  69  MPa,  £,  =  0.12%)  is 

=  58GPa(8.3Msi). 

At  this  stage  of  damage  development,  corresponding  to  the  proportional  limit  of  the 
stress-strain  response  of  the  0°  layer  (Fig.  2),  no  appreciable  damage  has  occurred  in  the 
0°  layer  (layer  1).  Simple  lamination  analysis  gives  the  following  axial  moduli  and  average 
axial  stresses  in  the  two  layers : 

El  =  127GPa(18.4Msi) 

E'2  =  22.8  GPa  (3.3  Msi) 

=  152  MPa  (22  ksi) 

<T2^-  27.6  MPa  (4.0  ksi), 

where  Ei  =  initial  modulus  of  0°  layer,  £2  =  degraded  modulus  of  90°  layer,  aij,,  (T2x  — 
average  axial  (x-direction)  stresses  in  layers  1  and  2,  respectively. 

Damage  initiation  and  development  in  0°  layer 

Prior  investigations  have  shown  that,  in  a  unidirectional  SiC/CAS  composite  under 
longitudinal  tension,  transverse  matrix  cracking  starts  at  an  applied  strain  of  approximately 
0.1%,  see  Daniel  et  al.  (1993).  This  value  is  slightly  below  the  proportional  limit  of  the 
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longitudinal  stress-strain  curve  of  Fig.  2.  Thus,  damage  in  this  layer  starts  near  the  satu¬ 
ration  point  of  cracking  in  the  90°  layer. 

The  initial  failure  in  the  0°  layer  consists  of  transverse  matrix  cracks  increasing  in 
density  with  applied  stress  up  to  a  minimum  crack  spacing  of  approx.  120  /rm  (0.0047  in.) 
at  a  stress  of  138  MPa  (20  ksi)  and  a  strain  of  0.37%,  as  shown  in  Fig.  6.  This  is  close  to 
the  minimum  crack  spacing  observed  before  in  a  unidirectional  specimen  under  longitudinal 
loading  and  corresponds  to  eight  fiber  diameters,  Daniel  et  al.  (1992). 

At  this  stage  of  damage  development  (3^=  138  MPa,  0.37%)  the  reduced  secant 
modulus  is 


E',=  37.3  GPa  (5.4  Msi). 


Assuming  no  increase  in  the  average  axial  stress  carried  by  the  damaged  90°  layer,  the 
following  axial  moduli  and  average  axial  stresses  are  obtained  in  each  of  the  two  layers: 


E\  =  97  GPa  (1 4.0  Msi) 
E\  =  7.5  GPa  (1.08  Msi) 
a^^  =  359  MPa  (52  ksi) 

=  27.6  MPa  (4.0  ksi). 


The  further  reduction  in  the  layer  2  modulus,  while  its  average  axial  stress  remains 
constant,  is  attributed  to  the  increased  opening  displacements  in  layer  2  macrocracks, 
following  crack  formation  and  multiplication  in  layer  1 . 

One  important  observation  is  that  the  in  situ  reduced  axial  modulus  E\  of  layer  1  is 
appreciably  higher  than  the  corresponding  modulus  of  the  unidirectional  material  at  the 
same  stage  of  damage  development,  i.e.  saturation  of  matrix  cracks  at  a  strain  of  s,  =  0.37% 
(see  Fig.  2).  This  could  be  attributed  to  the  fact  that  large  segments  of  layer  2  remain 
adhered  to  layer  1  and  bridge  the  matrix  cracks  in  layer  1,  thereby  increasing  its  effective 
in  situ  stiffness. 

At  the  same  point  of  damage  development,  there  is  evidence  that  extensive  fiber-matrix 
debonding  occurs  in  this  layer  near  and  after  matrix  crack  saturation,  see  Luo  et  al.  (1994), 
Wooh  and  Daniel  (1994).  As  in  the  case  of  the  unidirectional  material  under  longitudinal 
loading,  isolated  fiber  breaks  were  observed  before  saturation  of  transverse  matrix  cracking 
in  the  0°  layers. 


Damage  interaction  between  0°  and  90°  layers 

As  the  applied  stress  is  increased  above  138  MPa  (20  ksi),  the  matrix  cracks  in  the  0° 
layer  begin  to  propagate  into  the  90°  layer  and  connect  with  the  macrocracks  already 
developed  in  that  layer.  Because  of  the  difl'erence  in  crack  spacings  in  the  two  layers,  several 
cracks  in  the  0°  layer  connect  with  each  crack  in  the  90°  layer  forming  a  “delta-like”  pattern 
in  the  90°  layer  near  the  interface  (Fig.  7).  Each  macrocrack  in  the  90°  layer  corresponds 
to  five  to  six  transverse  cracks  in  the  0°  layer. 

With  increasing  load,  axial  (in  the  loading  direction)  cracks  are  generated  in  the  90° 
layer  near  the  interface  with  the  0°  layer  (Fig.  8).  These  cracks  (or  delaminations)  may  be 
due  to  high  interlaminar  shear  stresses  developed  at  the  tips  of  the  transverse  macrocracks 
in  layer  2,  and/or  high  interlaminar  normal  and  shear  stresses  existing  near  the  edges  of 
crossply  laminates  (edge  effect).  Additional  transverse  cracks  also  appear  in  layer  2. 

In  the  final  stages  of  damage  development  extensive  fiber-matrix  debonding  in  layer  1 
is  accompanied  by  fiber  fractures  and  sliding.  This  is  manifested  by  large  crack  openings  in 
layer  1  and  even  larger  (more  than  five  times  larger)  crack  openings  in  layer  2  (Fig.  9). 


SAS  32:3/4-F 


The  various  failure  mechanisms  and  the  entire  process  of  damage  evolution  are  cor¬ 
related  with  the  macroscopic  stress-strain  curve  of  the  [0/902]^  laminate  in  Fig.  10.  The 
effects  of  the  various  failure  mechanisms  on  stiffness  degradation  are  illustrated. 

Ultimate  failure  occurred  at  an  applied  stress  of  over  220  MPa  (32  ksi),  which  cor¬ 
responds  to  an  in  situ  longitudinal  strength  of  over  607  MPa  (88  ksi)  in  the  0  layer.  This 
is  much  higher  than  any  measured  longitudinal  strength  of  the  unidirectional  material  and 
clearly  shows  the  synergistic  effect  of  the  0“  and  90°  layers  of  the  laminate. 

4.  ANAl.YTICAL  PREDICTIONS 

The  formation  of  transverse  macrocracks  in  the  90°  layer  of  a  crossply  laminate  has 
been  discussed  by  several  investigators,  including  Reifsnider  et  al.  (1983),  Talreja  (1985), 
Hashin  (1985),  Ogin  et  al.  (1985),  Laws  and  Dvorak  (1988)  and  Lee  and  Daniel  (1990). 
The  analyses  above  were  aimed  primarily  at  polymer  matrix  composites,  but  since  they  are 
linear  they  can  be  applied  to  ceramic  matrix  composites  as  well.  The  analytical  model 
described  by  Lee  and  Daniel  (1990)  was  used  here  to  explain  the  observed  phenomena. 
This  model  was  shown  to  give  accurate  predictions  of  stiffness  degradation  and  crack 
spacing  in  the  90°  layer. 

Under  axial  loading  the  state  of  stress  in  each  layer  is  essentially  uniaxial.  The  first 
transverse  macrocracks  in  layer  2  occur  when  the  axial  (x-axis)  stress  in  that  layer  reaches 
the  strength  of  the  layer,  i.e.  the  transverse  tensile  strength,  Fj,,  of  the  unidirectional  material 
(Fig.  1 1).  As  discussed  before,  the  in  situ  transverse  strength  of  the  layer  is  somewhat  higher 
than  the  value  measured  by  testing  an  unsupported  unidirectional  layer.  The  applied 
laminate  stress  at  crack  initiation  in  layer  2  is 
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-ffL). 


(1) 


where 

Eh  \  E  h 

£  _  _J — ! - L_i  =  axial  modulus  of  undamaged  laminate, 

""  hf+hi 

El,  £2  =  longitudinal  and  transverse  moduli  of  undamaged  unidirectional 
material,  respectively, 

F2i  =  transverse  (to  the  fibers)  tensile  strength  of  unidirectional  material, 

ffjx  =  residual  axial  stress  in  90°  layer. 

The  above  residualstressisapproximatelyequal  to7  MPa  (I  ksi)  fora  1 100°C  (2000°F) 
temperature  difference  between  processing  and  room  temperature.  Substituting  the  known 
material  properties  and  the  value  £2,=  38  MPa  (5.5  ksi)  for  the  transverse  tensile  strength 
in  eqn  (I),  we  obtain  a  stress  of 


S,  =  28.8  MPa  (4.2  ksi), 

which  agrees  with  the  experimental  observation  of  initiation  of  transverse  macrocracks  in 
the  90°  layer  (Fig.  5). 

An  element  of  the  laminate  with  the  90°  layer  cracked  under  axial  loading  is  shown  in 
Fig.  11.  Of  particular  significance  to  damage  development  in  the  laminate  are  the  stress 
distributions  in  the  layers  and  the  stiffness  reductions.  The  stress  distributions  in  the  layers 
between  two  transverse  cracks  in  the  90°  layer  are  given  by  Lee  and  Daniel  (1990) 


ff,,.  = 


E. 


Fj/i,  cosh 


1  + 


“(I- 


EJi^  cosh  a 


^2 


(2) 


(3) 


(4) 


where  <r2;t=  average  (through  the  thickness)  axial  stresses  in  0°  and  90°  layers,  respec¬ 
tively,  Ti=  interlaminar  shear  stress  at  interface  between  0°  and  90°  layers, 

nhi  +  h2)E„  3G12G2J 

^~[hih2EiE2'  hiG22  +  h2Gi2\  ’ 

G,2,  G23=  in-plane  and  out-of-plane  shear  moduli  of  undamaged  unidirectional  material, 
=  axial  residual  stress  in  0°  layer  (calculated  as  — 13  MPa),  l2  =  crack  spacing  in  90 
layer. 
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This  crack  spacing  decreases  with  load  up  to  a  minimum  characteristic  value.  Exper¬ 
imentally,  a  minimum  crack  spacing  of  0.65  mm  (0.0255  in.)  was  observed  at  an  applied 
stress  of  d^  =  69  MPa  (10  ksi),  as  mentioned  before.  The  predicted  crack  spacing  at  this 
stress  can  be  obtained  by  setting  69  MPa  (10  ksi),  <T2x—  Pit—  38  MPa  (5.5  ksi)  and 
X  =  111!  (midpoint  between  cracks)  in  eqn  (3).  Then,  from  eqn  (3)  we  obtain 


where 


1 2  =  -cosh 


I- 


Plt 


El 


(6) 


For  the  numerical  values  above  we  obtain 

(/2)mm  =  0.68  mm  (0.027  in.). 

Expressions  for  the  reduced  axial  stiffness  E'2  of  the  90°  layer  and  5^  of  the  laminate 
are  given  as  a  function  of  crack  spacing  as  follows,  Lee  and  Daniel  (1990) : 


E'2^ 


(\-fS)(E2S.+  EM 


Eihi 

E,h] 


(1  -pWu 


and 


E'  =  E, 


1  + 


s,/i,  £|  J 


(7) 


(8) 


where 


At  the  applied  stress  level  of  69  MPa  (10  ksi),  where  transverse  crack  saturation  was 
noticed,  but  before  any  substantial  damage  occurred  in  the  0°  layer,  the  90°  layer  and 
laminate  moduli  are  calculated  as 


E'2  =  l9GPa(2.8Msi) 

£;=  55CjPa(8.0  Msi), 

which  are  in  good  agreement  with  the  experimental  values  of  E'2  —  22.8  GPa  (3.3  Msi) 
and  E'x  =  58  GPa  (8.3  Msi)  mentioned  before. 

The  average  stress  carried  by  the  0°  layer  then,  is 


146  MPa  (21.1  ksi). 


The  stress  necessary  to  initiate  transverse  matrix  cracking  in  the  0°  layer  is 
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li  =  6li 

Fig.  12.  Element  of  laminate  with  both  layers  cracked  up  to  their  saturation  levels. 


£.n 


(^int  ^\x 


where  =  matrix  modulus,  =  matrix  tensile  strength,  =  axial  micro-residual  stress 
in  matrix  of  O'"  layer. 

Neglecting  the  counteracting  residual  stresses,  tensile  micro-residual  stress  and 
compressive  macro-residual  stress  we  obtain  the  value 


154  MPa  (22.3  ksi). 

This  shows  that  crack  initiation  takes  place  in  the  0°  layer  at  a  slightly  higher  load  level 
than  that  at  which  crack  saturation  occurs  in  the  90°  layer.  This  also  corresponds  to  the 
observed  proportional  limit  of  the  longitudinal  stress-strain  curve  (Fig.  2). 

Cracking  continues  in  the  0°  layer  up  to  a  saturation  crack  density  where  (/i)min=  0.120 
mm  (0.0047  in.),  which  is  one-sixth  to  one-fifth  of  An  element  of  the  laminate  with 
both  layers  cracked  up  to  their  saturation  levels  is  shown  in  Fig.  12.  A  subelement  of  the 
90°  layer  at  the  interface  with  the  0°  layer  is  subjected  primarily  to  the  normal  axial  stress 
ffix  and  interlaminar  shear  stress  x,.  These  stresses  result  in  a  tensile  principal  stress  at  an 
angle  that  would  explain  the  formation  of  branch  cracks  in  the  observed  “delta*’  pattern. 


5.  SUMMARY  AND  CONCLUSIONS 

The  behavior  of  crossply  ceramic-matrix  composites  under  axial  tensile  loading  was 
studied.  The  material  investigated  was  SiC/CAS,  calcium  aluminosilicate  glass-ceramic 
reinforced  with  silicon  carbide  (Nicalon)  fibers.  Specimens  were  loaded  under  the  micro¬ 
scope  and  the  various  failure  mechanisms,  their  sequence,  interaction  and  overall  damage 
development  were  observed  in  real  time  and  recorded. 

The  first  stage  of  damage  development  consists  of  microcracks  in  the  90°  layer,  which 
then  develop  into  transverse  macrocracks.  These  macrocracks  increase  in  density  up  to  a 
saturation  limit,  with  a  minimum  crack  spacing  approximately  equal  to  the  90°  layer 
thickness.  Thereafter,  transverse  matrix  cracks  are  generated  in  the  0°  layer.  These  cracks 
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also  increase  in  density  up  to  a  minimum  crack  spacing  of  approximately  eight  fiber 
diameters.  In  the  third  stage  of  damage  development  there  is  an  interconnection  of  the  two 
sets  of  cracks.  Several  of  the  denser  0“  layer  cracks  connect  with  each  of  the  90°  layer  cracks 
in  a  delta-like  pattern.  This  is  finally  followed  by  delaminations  and  additional  cracking  in 
the  90°  layer  prior  to  ultimate  failure. 

Damage  development  was  correlated  with  the  macroscopic  stress-strain  response  of 
the  unidirectional  material  under  longitudinal  and  transverse  loading  and  of  the  crossply 
laminate.  It  was  noticed  that  first-ply  failure  in  the  90°  layer  occurs  at  a  higher  layer  stress 
than  the  transverse  tensile  strength  of  the  unidirectional  material.  This  may  be  attributed 
to  the  constraining  effects  of  the  0°  layer  which  tend  to  increase  the  in  situ  transverse  tensile 
strength  of  the  90°  layer. 

Following  crack  saturation  in  the  90°  layer,  cracking  initiated  and  increased  in  the  0° 
layers  at  approximately  the  same  strain  levels  where  similar  cracking  was  observed  in  the 
unidirectional  material  under  longitudinal  loading.  However,  the  0°  layer  stresses  at  the 
various  stages  of  damage  development,  were  higher  than  corresponding  stress  levels  in  the 
0°  unidirectional  specimen  for  the  same  damage.  This  was  explained  as  a  stiffening  or 
strengthening  effect  of  the  0°  layers  caused  by  the  attached  90°  layer.  Ultimate  failure  of 
the  laminate  clearly  showed  the  synergistic  effects  of  the  0°  and  90°  layers  of  the  laminate. 

Damage  development  in  the  90°  layer  was  compared  with  analytical  predictions.  A 
model  developed  previously  for  crossply  graphite/epoxy  composites  was  used  to  predict 
cracking  in  the  90°  layer  and  stiffness  degradation,  both  in  the  90°  layer  and  in  the  overall 
laminate.  The  predictions  were  in  very  good  agreement  with  experimental  observations. 

It  was  shown  that  the  overall  damage  development  is  characterized  by  different  scales 
at  different  stages.  The  relevant  scales  are  fiber  diameter  and  fiber  spacing  for  failure 
initiation  in  the  90°  layer;  layer  thickness  for  damage  saturation  in  the  90°  layer;  and  fiber 
diameter  for  matrix  cracking  in  the  0°  layer.  Finally,  there  is  a  scale  crossing  phenomenon 
in  the  interaction  between  the  damage  in  the  90°  and  0°  layers. 
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